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Vorwort 



Im Rahmen der Vorlesung „Einfiihrung in die Extremwerttheorie" von Prof. Dr. 
Martin Schlather lernte ich im Wintersemester 2008/2009 „ Brown- Resnick-Pro- 
zesse" als eine spezielle Klasse stationarer, max-stabiler stochastischer Prozesse 
kennen. 

Diese haben vielfaltige Anwendungen. So wird z. B. in [BdZ08j cin Brown-Resnick- 
Prozess (mit Standard-Frechet-Randern) zum Variogramm j((hx, /12)) = |(|^l| + 
I/12I) verwendet, um die raumliche Verteilung von extremen Regenfallen zu model- 
lieren. Aufgrund derartiger Anwendungen sind Brown-Resnick-Prozesse und ihre 
Simulation auch von starkem wirtschaftlichen Interesse, wie z. B. fur auf Monte- 
Carlo-Methoden basierenden Kalkulationen von Versicherern. 

Bei der Simulation von Brown-Resnick-Prozessen nach dem sich aus der Defi- 
nition kanonisch ergebenden Verfahren wachst jedoch der Approximationsfehler 
exponentiell mit der IntervallgroBe. Daher verhalt sich auch die Laufzeit entspre- 
chend schlecht, wenn der Fehler auf grofien Intervallen klein gehalten werden soli. 

Daraus ergeben sich die Problemstellung und das Ziel dieser Arbeit. Es sollen 
verschiedene Simulationsverfahren vorgestellt, untersucht und miteinander vergli- 
chen werden. Gesucht ist eine Anwort auf die Frage: Welches Verfahren ist - je 
nach Anforderungen an die Simulation - hinsichtlich des Fehlers „ optimal"? 

Die Arbeit gliedert sich dabei wie folgt: In Kapitel 1 werden die Grundbegriffe 
fur die Definition von Brown-Resnick-Prozessen bereitgestellt. Diese erfolgt in 
Kapitel 2 ebenso wie die Vorstellung stochastisch aquivalenter Darstellungen der 
Prozesse. Das Kapitel basiert im Wesentlichen auf dem Artikel „ Stationary Max- 
Stable Fields Associated to Negative Definite Functions" von Sachar Kablutschko, 
Martin Schlather und Laurens de Haan (im Weiteren als |KSd09| bezeichnet) sowie 
weiteren Ideen der Autoren. 

Anhand der unterschiedlichen Darstellungen ergeben sich verschiedene Simulati- 
onsverfahren fur Brown-Resnick-Prozesse. Diese werden in Kapitel 3 naher be- 
schrieben und der bei den Verfahren auftretende Approximationsfehler wird ab- 
geschatzt. SchlieBlich werden die Ergebnisse der Simulationsverfahren in R in Ka- 
pitel 4 dargestellt und verglichen. Programmcode, Dokumentation und eine detail- 
lierte Auflistung der Ergebnisse der zugrundeliegenden Simulationsstudie finden 
sich im Anhang. 

Sehr herzlich danken mochte ich Prof. Martin Schlather fur meine Ausbildung in 
vielfaltigen Bereichen der Stochastik, insbesondere stochastischen Prozessen und 
Extremwerttheorie, fur die Vermittlung wesentlicher Programmierkenntnisse in R, 
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die Bereitstellung des Packages RandomFields sowie fur die spannende Fragestel- 
lung und die hervorragende Betreuung der Arbeit. 

Weiterhin gilt mein Dank Dr. Sachar Kablutschko fur seine Hilfe bei Messbar- 
keitsproblemen sowie seine Bereitschaft zur Zweitkorrektur der Arbeit. 

Gottingen, im Juli 2009 
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= Gleichheit in Verteilung 

Af([jL, S) Normalverteilung mit Erwartungswert \x und Kovarianzmatrix S 

Exp(A) Exponentialverteilung mit Parameter A 

U(M) Gleichverteilung auf einer Menge M 
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1 Einleitung 



1.1 Zufallsfelder 

In diesem Kapitel sollen zunachst die Begriffe, die fiir die Definition von Brown- 
Resnick-Prozessen benotigt werden, erlautert werden. Dabei handelt es sich u. a. 
urn grundlegende mogliche Eigenschaften von Zufallsfeldern. 

Definition 1.1. Sei T eine Indexmenge. Eine Familie von reellwertigen Zufalls- 
variablen {Z(t), t G T} heifit stochastischer Prozess. Im Fall T = R rf nennt man 
{Z(t), t G T} auch d-dimensionales Zufallsfeld. 

Definition 1.2. Ein Zufallsfeld heifit Gaufisch genau dann, wenn alle endlich- 
dimensionalen Randverteilungen multivariat normal sind. 

Definition 1.3. Ein Zufallsfeld Z(-) auf ~R d heifit 

1. stationdr (im engeren Sinne) genau dann, wenn alle endlich-dimensionalen 
Randverteilungen translationsinvariant sind, d. h. 

P(Z(xi +x) G Bi, . . . , Z(x n + x)e B n ) = P(Z(xi) G B\, . . . , Z(x n ) G B n ) 
fiir alle x, Xi G R d , Bi G B d , i G {1, . . . , n}, n G N. 

2. stationdr (im weiteren Sinne) oder schwach stationdr genau dann, wenn 
Erwartungswert KZ(t) und Kovarianz Cov(Z(x + t),Z(t)) fiir alle x G M d 
nicht von t G M d abhangen. 

Die Funktion C : R d -)• R, 

C(/i) = Cov(Z(/i),Z(0)) 

heifit Kovarianzfunktion. 

Bemerkung 1.4. Da eine multivariate Normalverteilung eindeutig durch Er- 
wartungswert und Kovarianzmatrix bestimmt ist, sind in dem Fall, dass Z{-) ein 
Gaufisches Zufallsfeld ist, die Begriffe „stationax im engeren Sinne" und „stationar 
im weiteren Sinne" aquivalent. 

Definition 1.5. Ein Zufallsfeld Z(-) auf R rf heifit intrinsisch stationdr genau 
dann, wenn fiir jedes h G R d das Zufallsfeld {Z(x + h) - Z(x), x G R d } schwach 
stationar ist. 

Die Funktion 7 : R d -)■ [0, 00), 

7 (fc) = iE(Z(fc) - Z(0)) 2 

heifit ( Semi- ) Variogramm. 
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1. Einleitung 



Bemerkung 1.6. Aus der Stetigkeit des Variogramms im Ursprung folgt die 
Stetigkeit des Variogramms auf R , derm es gilt fur beliebiges x G R d : 

lim E(Z(x + h)- Z(0)) 2 - E(Z(x) - Z(0)) 2 

= lim E(Z(x + h)- Z(x)f + 2E(Z(x + h) - Z(x))(Z(x) - Z(0)) 

+ E(Z(x) - Z(0)) 2 - E(Z(x) - Z(0)) 2 
= lim E(Z(h) - Z(0)) 2 + 2E(Z{x + h) - Z(x))(Z(x) - Z(0)) 

h— >oo 

< lim E(Z(/i) - Z(0)) 2 + 2y / E(Z(x + h) - Z(x)) 2 ■ E(Z(x) - Z(0)) 2 

h— >oo 

nach der Cauchy-Schwarz-Ungleichung 

= 0, falls lim E(Z(h) - Z(0)) 2 = 0. 
h— >0 

Somit ist ein Variogramm 7 genau dann stetig auf M. d , wenn es im Ursprung stetig 
ist. 

Definition 1.7. Sei {Z(t), t G R d } ein Zufallsfeld auf R d mit Standard-Gumbel- 
Randern, d. h. es ist P(Z(i) < x) = exp(— e~ x ) fiir alle t € M. d , x G R, und seien 
Z(-) Z(-), i G N. 

Dann heifit Z(-) max-stabil genau dann, wenn fiir alle n G N die Beziehung 

max Zj(-) — logn = Z(-) 
i=l, ...,n 

gilt. 



1.2 Poisson-Punktprozesse 

Weiterhin treten im Zusammenhang mit Brown-Resnick-Prozessen sogenannte 
„Poisson-Punktprozesse" auf, die auf „Polnischen Raumen" definiert sind. 

Definition 1.8. Sei (fl, m) ein metrischer Raum und T die von m induzierte 
Topologie. Dann heifit f2 Polnischer Raum genau dann, wenn beziiglich m 
vollstandig ist und T eine abzahlbare Basis besitzt. 

Die von der Topologie T erzeugte a- Algebra, <t(7"), heifit Borel-a -Algebra. 

Definition 1.9. Sei S ein lokal kompakter Polnischer Raum mit Borel-o"- Algebra 
S, versehen mit einem lokal endlichen Mafi [i. Ein zufalliges Zahlmafi IT auf (S,S) 
heifit Poisson-Punktprozess mit Intensitatsmafi /1 genau dann, wenn 

1. IT(j4i), . . . , U(A n ) sind stochastisch unabhangig fiir alle A±, . . . , A n G S mit 
Ai n Aj = fur i / j, i,j G {1, . . . , n} und n G N, 

2. I1(j4) ist Poisson-verteilt mit Parameter fi(A) fiir alle A G 5. 
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1.2. Poisson-Punktprozesse 



Das nachfolgende Lemma behandelt die „Erweiterung" von bestimmten Poisson- 
Punktprozessen um Zufallsvariablen. Eine solche wird im weiteren Verlauf der 
Arbeit noch mehrfach auftreten. 

Lemma 1.10. Sei X^ieN^i P°^ sson -P un ktprozess auf R mit Intensitdtsmafi 
Xe~ x dx, A > 0. Seien weiterhin Y±,Y2,... unabhdngig identisch verteilte Zufalls- 
grofien (unabhdngig von {Xi}i & fq) mit Werten in einem Polnischen Raum S. Sei 
P das zu Yi gehdrige Wahrscheinlichkeitsmafi. 

Dann ist LT := X^ieN $(Xi,Yi) e ^ n Poisson-Punktprozess auf KxS mit Intensitats- 
mafi \e~ x dxP(dy). 



Beweis. Zunachst bemerken wir, dass das MaB Ae x dxP(dy) lokal endlich ist, 
denn fur jede beschrankte Borel-messbare Menge A x B C R x S ist 

[ \e- x dx¥(dy) < [ X e - X dx¥(dy) = [ \e~ x dx < oo, 
JAxB JAxS Ja 

da A beschrankt ist. 

Seien nun [a, b] C R, A C S Borelmengen. Dann gilt fur k G N: 
P(n([a,6] x A) = k) 

oo 

= ]T P(n([a, b] = n)¥(#{i €{l,...,n}:Yi€A} = k) 

n=k 

= Y exp(-A(e" a - e - fe ))l(A(e- a - e~ b )) n ( "V(A) fe (l - ¥(A)) n ~ k 

= exp(-A(e- a - e~ b )) (p(,4)A(e- a - e^ 6 )) * 

■ f; We-"-e-^ (1 _ p(A)r 

n=0 



i exp (-P(^)A(e" a - e~ 6 )) (p(^)A(e" a - e" 
^ exp ( - J A y e _ *dzdP(i/) J IJ X f e^ x dxdF(y) J 



,,6] / J[a,6] y 

Die Anzahl von Punkten in einer Borel-Menge ist also Poisson-verteilt mit der 
behaupteten Intensitat. 

Die Unabhangigkeit der Anzahlen von Punkten in Mengen, die in der ersten Kom- 
ponente disjunkt sind, folgt bereits aus der Tatsache, dass J2ieN em Poisson- 
Punktprozess ist und dass Y±,Y2, . . . unabhangig sind. 

Es geniigt daher, dies fur Mengen der Form [a, b] x A\, [a, b] x A2, ■ ■ ■ , [a, b] x A n 
zu zeigen mit disjunkten Borel-Mengen Ai, A2, ■ ■ ■ , A n C S und [a, b] C R. Seien 
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1. Einleitung 



dazu fei, . . . , k n G N, k\ + . . . + k n = k und = A(e a — e ). Dann gilt 
P(fi(Ai x [a,6]) = fci,...,n(A n x [a, 6]) = fc n ) 

j = fc J ' U \ i = l / 8=1 *" 

n 

= II rT ( P ( A 0Aa, b )" 1 exp (-A a , 6 P(^)) 

1=1 

n 

= [J p (n(A x [a, 6]) = A^), 

i=l 

also die zu zeigende stochastische Unabhangigkeit. 

Somit ist II ein Poisson-Punktprozess mit dem angegebenen Intensitatsmafi. □ 

Ebenfalls werden wir den folgenden Transformationssatz fiir Poisson-Punktpro- 
zesse benotigen. 

Satz 1.11. Seien S, T Polnische Raume mit Borel-a-Algebren S bzw. T und 
II = X^ieN e ^ n Poisson-Punktprozess auf (S, S) mit lokal endlichem Inten- 
sitatsmafi [i und f : S — > T messbar so, dass auch das Bildmafi fi* = fio/ _1 lokal 
endlich ist. Dann ist durch X^ieN^/(^i) e ^ n Poisson-Punktprozess auf (T,T) mit 
Intensitatsmafi [i* gegeben. 

Beweis. siehe |Kle06j . Satz 24.16. □ 
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2 Brown-Resnick-Prozesse 



Bevor Brown-Resnick-Prozesse definiert werden konnen, sind zunachst noch ei- 
nige Vorbemerkungen notwendig: In diesem Kapitel werden wir u. a. Poisson- 
Punktprozesse auf dem Raum C(]R d ), also dem Raum der reellwertigen, stetigen 
Funktionen auf betrachten. Dazu versehen wir C(R d ) mit der a- Algebra C, 
die von den Mengen der Form 

C tl ,...,t m (B) := {/ G C(R d ) : (f(h), f(t m )) G B}, 

mit tx, ■ ■ ■ , t m G M. d , m G N und B G B m erzeugt wird. Diese a- Algebra ermoglicht 
es uns, Aussagen iiber Funktionswerte an endlich vielen Stellen zu treffen, was sich 
als sehr niitzlich herausstellen wird. 

Nun ist Qk := C([— k, k] d ), versehen mit der Metrik irik gegeben durch 

m k {f,g)-= sup \f(x)-g(x)\, 

x£[-k,k] d 

ein Polnischer Raum. Dabei ist durch 

T k :={Bf\f): e G Q>o, / G IIq}, 

wobei IIq die Menge aller Polynome auf W 1 mit rationalen Koeffizienten und 

Be (/) := {g G '■ mk{f,g) < e} sei, nach dem Satz von Stone- WeierstraB 
(reelle Fassung) eine abzahlbare topologische Basis gegeben. 

Damit ist auch der Produktraum £1 := X'kLx £lk m it der Metrik 

oo 

tx 1 + m * 

ein Polnischer Raum. 



Lemma 2.1. Die von m auf £1 induzierte Topologie ist die Produkttopologie aus 
den von den Metriken aufVL^, k G N, induzierten Topologien. 

Beweis. Es geniigt zu zeigen, dass die Basis der einen Topologie jeweils in der 
anderen Topologie enthalten ist. 

Eine Basis der Produkttopologie bilden die Mengen der Form X keN B ( S{f k ) mit 
fk G Slfc, wobei nur endlich viele e k endlich seien und wir B^\f) := definieren. 
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2. Brown- Resnick-Prozesse 



Fur g £ X ke ® B^ (f k ) mit g = (51,52, •••) gibt es nun 6 = (61,62,...) mit 
< 6k < 00, falls e k < 00, und 6 k = 00, falls e k = 00, so, dass 

X B^(g k ) C X B^(f k ). 
fceN ken 

Sei e := min {feeN . Sk<oo} (2 _fe_1 -min^l)), dann folgt - wegen ± x+mjA) < 
m(-, •) - fiir h = (hi, hi, ■ ■ ■) £ -Be(fl') := {/i £ ^ : m(/t, 5) < e}, dass 

m(g,h)2 k m(g, h)2 k 

m k (g k ,h k ) < - — — -r < j < 6 k , 

l-m(g,h)2 k ± 



falls 6 k < 00. Somit ist 



s e ( fl ) c X c X B^(f k ), 

k&N feeN 



d. h. X Be$(fk) ist auch offen beztiglich der von m erzeugten Topologie. 

Sei andererseits nun ein Basiselement B £ (f) der von m induzierten Topologie 
gegeben. Dann existiert fiir jedes g £ B e (f) ein 6 > mit B$(g) C B e (f). Wahle 
nun iV £ N so grofi, dass 2 _Ar < §. 

Dann gilt fiir jedes /i £ £$,(#1) x Bj^(ga) x . . . x B^(g N ) x x£ljv +1 fi fc , dass 

JV 00 iV r XX 

m( 5 ,/ i )<^2- fe m fe ( 5fe) ^)+ 2 2^<^2- fc - + 2- JV <- + - = 5. 

fe=l j=7V+l fc=l 

Also ist Bff 2 ( gi ) x B$Gfe) x ... x flgjM x Xf =N+1 n k c Ba( ff ) C B E (f), d. h. 
jedes Element g £ B e (f) hat eine Umgebung aus der Basis der Produkttopologie, 
die in B e (f) enthalten ist. Somit lasst sich B e (f) als Vereinigung von Basisele- 
menten der Produkttopologie schreiben und ist daher in dieser enthalten. □ 

Da wir C(]R rf ) mit der abgeschlossenen Teilmenge 

F :={(h, f 2 , ...) £ n : fn+i\ { _ nMd = fn Vn £ N} C fi 

identifizieren konnen, ist auch C(]R d ) mit der Metrik m ein Polnischer Raum. 
Nun ist der Zusammenhang zwischen der Borel-cr-Algebra von C(M d ) und der 
er-Algebra C zu klaren. 



Proposition 2.2. Sei r die von m auf C(R ) induzierte Topolgie. Dann gilt 
a(r)=C. 

Beweis. Wir stellen fest, dass C bereits von den Mengen der Form CVt,...,t m (B), 
wobei t\,...,t m £ W 1 , m £ N und B C M m offen seien, erzeugt wird. Da B offen 
ist, existiert zu jedem / £ Ct lr .. ) t m (B) ein e > mit 

C tl ,.. . ((/(ti) - e, f(h) + e) x ... x (/(i m ) - e, f(t m ) + £ )) c C tl ,...,t m (B) i 
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denn diese Teilmenge ist ein offener Ball um / in der Supremumsnorm auf R m . 

Wahlen wir nun N G N derart, dass t\, . . . , t m G [-N, N] d ist, so gilt nach Kon- 
struktion 

([XB^g)x x n k ) HF cC tl _ tm (B). 

\ \k=l k=N+l J j 

Also gibt es zu jedem Element aus Ct 1: ... : t m (B) eine in Cti,...,t m (-B) enthaltene offe- 
ne Umgebung aus der Basis der Spur der Produkttopologie auf F. Diese Topologie 
ist nach Lemma HO gerade r. Also ist jede erzeugende Menge Ct 1 ,....t m {B) G r und 
damit gilt auch C C o"(r). 

Da 

T := | ( X flJH/*)] HF: e fc G Q> U {oo}, / fc G IIq, 
[ \fceN / 

wobei nur endlich viele e k endlich seien > 



nach obigen Ausfiihrungen und Lemma 12.11 eine abzahlbare Basis von r ist, gilt 
a"(J-") = ct(t) und es geniigt G G C fiir alle G £ Jzu zeigen. 
Sei nun q^jq^, ■ ■ ■ eine Abzahlung von <Q> d n [— A;, A;] a! . Dann ist 

( X BW(f k )) nF 

5 G C(R d ) | VA: G N 3n k G N Vj G N : | 5 (g. (fc) ) - f k (qf ] )\ < e k - — 

= n u nc,.((A(«f)-« + ^.A(#») + «-^))€ C 

- wie behauptet. □ 

Somit ist C(]R d ) ein Polnischer Raum mit Borel-c-Algebra C. Da C(M. d ) zudem 
lokal kompakt ist, konnen wir im Folgenden auf (C(M. d ),C) Poisson-Punktprozesse 
definieren. 
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2. Brown- Resnick-Prozesse 



2.1 Brown-Resnick-Prozesse als max-stabile Prozesse 

Fur den Beweis des Satzes 12.41 der gleichzeitig den Begriff Brown-Resnick-Prozess 
definiert, benotigen wir das folgende Lemma. Dieses lasst sich leicht nachrechnen, 
weshalb hier auf seinen Beweis verzichtet werden soil. 



Lemma 2.3. Seien 

(Xi,X 2 ) ~ TV 



fJ-i \ f E 12 



Dann gilt: 

X 2 \(X 1 = xi) ~ A/"(/i 2 + S 2 iS5" 1 1 (xi - mi),S 2 2 - S2iS5" 1 1 Si 2 ). 



Der nachfolgende Satz geht auf Brown und Resnick zuriick, die ihn fur die Brown- 
sche Bewegung als zugrundeliegendem Prozess (d. h. fur den Spezialfall 7(^1) = 
\\h\) bewiesen (vgl. [BR77]). 

Satz und Definition 2.4. Sei W(-) ein pfadstetiges, intrinsisch stationares 
Gaufisches Zufallsfeld auf M rf mit Variogramm j, Erwartungswert KW(t) = 
und Varianz a 2 (t) := Var(W(t)), t E R d . 

Weiterhin seien Wi(-) ~ W(-), i E N unabhdngig identisch verteilt, sowie 

IT = Yl ?>Xi e * n t ' on diesen stochastisch unabhdngiger Poisson-Punktprozess auf M 

mii Intensitatsmaft e~ x dx. Dann ist durch 

Z(t) := max (Xi + Wdt) teR d 
ieN \ 2 / 

ein stationdrer max-stabiler Prozess auf R d mit Standard- Gumbel-Randern gege- 

ben. Die Verteilung von Z hdngt nur von 7 a&. 

Z(-) /teijSi Brown-Resnick-Prozess zum Variogramm 7. 

Beweis. Seien ti, . . . , t m E M d , yi, . . . y m E M, m E N beliebig und P 4 ~ f ~ sei 

das zum Zufallsvektor (W(ti), . . . , W(£&)) gehorende Wahrscheinlichkeitsmafi. 
Wir definieren nun fur i = 1, . . . ,m die GroBen of := cr 2 (ij) sowie die Zufallsva- 

riablen & := J^(tj) - ^- Dann ist & ~ Af(-4, of)- 



Sei 



r y := Covfe - - 6) = E((W(*i) - W(ti)) • (Wfo) - W(*i))) 
2 



,-E{((W(ti) - W(h)) - (W{t 3 ) - W(h))) 2 



- {Wit,) - W{ tl )f - {W{t 3 ) - W(h)) 2 } 
l(U - tj) + 7(*i - h) + j(tj - h) (2.1) 
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2.1. Brown- Resnick-Prozesse als max-stabile Prozesse 



und T := (Tij). 



l 3 H=2,...,m, j=2,...,m' 



Weiterhin seien /x 2m := 0E(Cj - Ci))j=2,...,m = 
und g := (Cov(£j - fi,£i))j=i,.„, m . 
Nun ist 

/ e 2 -a \ /-ii o 



Cm - Cl 



-1 

V i o 



j=2,...,m 



6 



o 1 
/ 



V Cm / 



und (Ci, • • • , Cm) ist normal verteilt. 



Daher ist auch (C2 — Ci> • • • i Cm — Ci> Ci) normalverteilt und nach obigen Uberle- 
gungen gilt 



/ C2-C1 \ 

Cm - Cl 

V Ci 



Dann ist nach Lemma EI 



r 9 



1 ' \ g T a\ 



Cl I (Ci-Cl)i=2,...,m ~ M (-^ +/r- 1 ((Ci - Cl)i=2,...,m " M2m),*? " 

(2.2) 

Weiterhin ist $ := SiLi ^(Xi,Wi(-)) nac h Lemma 11.101 ein Poisson-Punktprozess 
auf (M. x C(R d ),B x C) mit Intensitatsmafi e~ x dx x P und unter Verwendung der 
Aquivalenz 



Z(h) <yx,.. .,Z(t m ) < y Tl 



Xi < min y k - Wi(t k ) + 

k=l,...,m \ Z 



Vz € N 



$ N(s,iu) :x> k mm m (y k - w(t k ) + a 2 {t k )/2) X) 



folgt daraus 



log(P(Z(ti) <yi,...,Z(t m ) <y m )) 

e~ x dx f tl ,...,t m (d(wi, . . .,w m )) 

minfe=i,..., m (yfe-Wfe+o- 2 (t fc )/2) 

max e-^-^+ CT2 ^)/ 2 )p tl ,... )tro (d( U ; 1 , . . .,w m )) 

k=l,...,m 
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2. Brown- Resnick-Prozesse 



E| max e " fefc " ?fe) J = E ( e ?1 max e ^ k ~^~ yk 

,k=l m J \ fc=l,...,m 



= E ^E ((& t max m e^-^-^ | - £i)i=2,...,m)) 
= E ( max e & - ?l) - yfc E ( e ?1 I (£,- - £1)^ j 
^Ef max e^-^-exptg^-^J-^j-^-lg^g)) 

\k=l,...,m \ Z I J 

da U ~ A/"0, (T 2 ) E(e c/ ) = e^+^r 
= / max {e-f^e^-^j-expf/r- 1 ^-^)- Vr- 1 ^ 

jRm-i k=2,...,m \ I J 



(2 7 r) I Vdet(r) V 2 

/ max { e -^,e^-^} — \ 

J Rm -i k=2,...,m (2tt)^~ 



exp ( --(w - /i2m) T r 1 (u;-/i2m)) dw 2 ■ ■ ■ dw r , 



(2ir)— det(r) 

-xC 1 " 7 _ M2m - 5) T r _1 (w - //2m ~ 5) ) 



Mit 

(//2m + = - + E((^ - E£j - 6 + - Eft)) 

= - ^E (-(ft - Eft) 2 + (0 - Eft) 2 - 2(ft - Eft)(ft - Eft) + 2(6 - Eft) 2 ) 
= - ^E(ft ~Mj -6 + Eft) 2 = -^W,) - W^)) 2 = -7(tj - h) 

erhalt man also 



-log(¥(Z(t 1 )<y 1 ,...,Z(t m )<y m )) 

[ max { e -y\e Wk - yk }-—— ^ 

J Rm -i k=2,...,m l ' (2vr)( m ^ 1 )- 



!)/ 2 det(r) 

exp ( --(it; + 7(t. - ti)) T T~ l (w - j(t. - ti)) ) dw 2 ■ ■ ■ dw Ti 



In dieser Gleichung treten nach (|2.ip niir Term© dcr Form ^(tj — ^)?j=i,...,m 

auf, 

d. h. Z{-) ist stationar und die Verteilung hangt nur von 7 ab. 
Im Fall m = 1 ist das Integral degeneriert und wir erhalten 

-log(P(Z(t) <y) = e-w, 

d. h. Z(-) hat Standard-Gumbel-Rander. 
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2.1. Brown- Resnick-Prozesse als max-stabile Prozesse 



beien nun 

^u.i.v. Z(-) fu.r I G {1, . . . ,n\. Dann gilt 
- log(P( max Z^ih) - logn < Vl , . . . , max Z®(t m ) - logn < y m )) 

1=1, ...,n l=l,...,n 

= - log(P(Z(ti) < yi + log n, . . . , Z(i m ) < y m + log n)) • n 

1 



/ max { e -fi- 1 °g™ )e «'fc-^- lo s™} 

Jjm-i k=2,...,m 



(2^)( m " 1 )/ 2 det(r) 



-^(-w + 7(t. - ti)) T r 1 (w - 7(t. - ii)) 



- / max {e^'.e^^} 

77, Jjgm-i k=2,...,m 



1 



( 2vr )(m-l)/2 det ( r ) 



• exp ^-^(io + 7(t. - ti)) T r 7(t. - ti))J dw 2 ■ ■ ■ dw m ■ n 

= - log(P(Z(ti) < yi, . . . , Z(t m ) < y m )), 

also max; = i n - wie behauptet. □ 

Bemerkung 2.5. 

1. Im Folgenden wird mit Z{-) stets der oben definierte Brown-Resnick-Prozess 
bezeichnet. 

2. Sei W(-) ein intrinsisch-stationares Gaufisches Zufallsfeld auf M. d mit Erwar- 
tungswert und Variogramm 7. Existieren a, f3,Ck > der art, dass 

E||W(fc) - W(0)\\ a < C k \\h\\ d+ P fiir alle h G fc G N, (2.3) 

so gibt es nach dem Satz von Kolmogorov-Chentsov (Verallgemeinerung auf 
R d , vgl. Abschnitt I.I.I. 3 in |BS02j ) eine pfadstetige Modifikation von W{-). 
Ob diese Voraussetzung erfiillt ist, hangt nur vom Variogramm 7 ab, denn es 
ist W(h) — W(0) ~ M(0,2j(h)). Wir konnen daher fiir die entsprechenden 
Variogrammtypen die Pfadstetigkeit ohne Einschrankung voraussetzen (vgl. 
hierzu Abschnitt 2.3.1 in |CD99| V 

3. Definition und Satz 12.41 lassen sich auch auf nicht-pfadstetige Zufallsfelder 
erweitern. Aus technischen Griinden wird in dieser Arbeit jedoch darauf 
verzichtet. 

4. Wir konnen auch Brown-Resnick-Prozesse mit Standard-Frechet-Randern 
(d. h. F(Z(t) <x)= exp(-i) Vx G R, t G R d ) erhalten durch 

Z(t) = max exp ( X { + WAt) - ^-^] , t G R d . 
ieN y 2 J 

Ebenso erreichen wir Standard-Frechet-Rander in den weiteren Satzen, in- 
dem wir auf die dortigen Prozesse vor der Maximumsbildung die Exponen- 
tialfunktion anwenden. 
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2. Brown- Resnick-Prozesse 



2.2 Stochastisch aquivalente Darstellungen mit zufalligen 
Translationen 



Satz 2.6. Seien W wie in Satz \2.4\ W% W, i G N. Sez Q ein Wahrschein- 

lichkeitsmaj] auf R d und II = X^p£"i,-Hi) von ^ G N} unabhdngiger 

Poisson-Punktprozess auf KxK d mrf Intensitatsmafi e~ x dx Q(dh). Dann ist 

Z(t) := max (Xi + Wi(t- Hi) i G R d 
ieN \ 2 J 

ein Brown- Resnick-Prozess zum Variogramm 7, d.h. Z = Z . 

Beweis. Seien t±, t m G R d , y\, . . . y m G R, m G N beliebig und Pr ,~ 
sei das zum Zufallsvektor (W(ti), . . . , W(t k )) gehorende Wahrscheinlichkeitsmafi. 
Nun ist := X^ieN $(Xi,Hi,Wi) n ach Lemma [1.101 ein Poisson-Punktprozess auf 
(R x R d x C{R d ), B xB m x C) mit Intensitatsmafi e" x <ix Q(d/i)P(dW) und es 
gelten die Aquivalenzen 

Z(t\) < yx, . . . , Z(t m ) < y m 

X t < min (y k - Wi{t k - Hi) + a2 ^~ H i) \ \/z G N 

fc=l,...,m y A J 

<^ §(Ux,h,w) G R x R d x C(R d ) : 

x > min (y k - w{t k - h) + - ^ )) = 0. 

fe=l,...,m y A J ) J 

Damit ergibt sich 

-log(F(Z(t l )<y 1 ,...,Z(t m )<y m )) 

/ , tr 2 (t k -h) \ 
muifc=l,...,n»l Kfc-tWfcH § I 

Pt 1 -h,... ) t m -h(d(wi, ■ ■ • , »m)) Q(d/i) 
= - log ( f F(Z(h -h)< yi ,... Z(t m -h)< y m ) Q(dh) 

= - log f / P(Z(*i) <i/i,... Z(t m ) < y m ) Q(d/i) 

da Z(-) stationar ist nach Satz [2 
= - log(P(Z(ti) <!/!,... Z(t m ) < y m )), 

da Q e i n Wahrscheinlichkeitsmafi ist. 
Also gilt Z = Z - wie behauptet. □ 

Dieser Satz ermoglicht uns stochastisch aquivalente Darstellungen des Brown- 
Resnick-Prozesses. Zwei Spezialfalle werden in den folgenden Korollaren betrach- 
tet. 
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2.3. Brown- Rcsnick-Stationaritat 



Korollar 2.7. Seien W wie in Satz\2J\ W\ 3) ~ u .i.v. W, i 6 N, j = 1, ...,n. 

Weiterhin seien IIW = „yj, j = 1, . . . ,n stochastisch unabhangige Poisson- 

i 

Punktprozesse auf R (unabhdngig von {W i (•), i 6 N, j € {1, . . . , n}}) mit 
Intensitdtsmaji — ■ e~ x dx sowie hi, . . . , /i n G R rf /est. Dann ist 

Zi (t) := max max f X (i) + (t - fc f ) - ^ , (£K d 

j=l,...,n ieN \ 4 4 J 2 / 

em Brown- Resnick-Prozess zum Variogramm ^ , d. h. Z\ = Z . 

Beweis. Nach Konstruktion ist X^eN fis ) ^ r J e< ^ es J = 1 , - - - , ^ ein Poisson- 

Punktprozess auf R x R rf mit IntensitatsmaB ^ e~ x dx 8^. (dh). Diese Prozesse sind 
fur j G {1, . . . , n} stochastisch unabhangig. 

Daher ist auch X)j=i SieN ^fx^ 5 h ) e * n P°i sson -P un ktprozess und das zugehorige 
IntensitatsmaB ist ^e~ x dx(J2 7 j=i$h j )(dh). Mit Q = ^X)?=i^/y daher die 

Aussage des Korollars aus Satz 12.61 denn Q ist ein WahrscheinlichkeitsmaB auf 
R d . □ 



Korollar 2.8. Seien W, Wi, i G N wie in Satz \2.4\ /cK d ein endliches Intervall. 
Weiterhin sei II = z^^(Xi,Hi) davon unabhdngiger Poisson-Punktprozess auf 
1x1 mit Intensitdtsmaji \I\ e~ x dx dh. Dann ist 

Z 2 {t) := max f X< + W t (t - Hi) - ° 2 ( t - H i) \ t G R d 

ein Brown-Resnick-Prozess zum Variogramm 7, d. h. Z2 = Z. 

Beweis. Es ist durch Q(dh) := rkdh das zur Gleichverteilung auf I gehorende 
WahrscheinlichkeitsmaB gegeben. Damit folgt die Aussage aus Satz 12.61 □ 



2.3 Brown- Resnick-Stationaritat 

Der folgende Abschnitt liber den Begriff der „Brown-Resnick-Stationaritat" ori- 
entiert sich an Kapitel 2 aus |KSd 09j . 

Sei {£(t), t G R rf } ein stochastischer Prozess mit der Eigenschaft 

E(e m ) < 00 Vt G R d (2.4) 

und seien ~ u .i.v. £(")> i G N. Sei X^eN^Q e ^ n davon unabhangiger Poisson- 
Punktprozess auf R mit IntensitatsmaB e~ x dx. 
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2. Brown- Resnick-Prozesse 



Definition 2.9. £(•) heifit Brown- Resnick- stationar genau dann, wenn der Pro- 
zess {r)(t), t £ W 1 }, gegeben durch 

T)(t) :=max(Xi+ti(t)), 

stationar ist. 



Bemerkung 2.10. Sei W(-) ein pfadstetiges, intrinsisch stationares Gaufisches 
Zufallsfeld auf M. d mit Variogramm 7, Erwartungswert und Varianz <J 2 (t). 
Wahlen wir nun f (i) = W{t) - , dann ist 

Ee «*)= /■ e «-4 a _i e -K^)) 2 ^= /" — i e'K 2 ^) d«; = l 

7r v^7r<r(t) y/2Ku{t) 

fur alle i £ Damit ist £(•) Brown-Resnick-stationar, denn es ist fy(-) = Z'(-) 
stationar nach Satz 12.41 



Wir wollen nun eine aquivalente Definition der Brown-Resnick-Stationaritat ge- 
ben. Dazu betrachten wir den Raum C(M. d ) versehen mit der Borel-a- Algebra C. 
Nun konnen wir nach Satz ll.lTI als Poisson-Punktprozess auf C(M d ) 

auffassen und das zugehorige IntensitatsmaB A ist das Bildmafi von e~ x dx x P 
(wobei P das zur Verteilung von £(•) gehorende Wahrscheinlichkeitsmafi sei) unter 
der messbaren Abbildung 

vr : M x C(R d ) -> C(R d ), (x, £(•)) x + £(•)» 

also 

A(C tl! ... )tm (B))= /" e-* f F tl! ... itm (d Wl ,...dw m )dx 

JR JB-(x,...,x) 

= ^B-(x,...,x)(wi,---,w m ) e~ x dxP tu ... !tn (dwi,...,dw m ). 

JR d JR 



Fiir t £ R d , z £ Z betrachten wir die Menge A t>z := {/ £ C(R d ) : f(t) > z}. Dann 
ist 

A(A t>z )= [ e- x nm>z-x)dx v '- = = x [ e y- z nm > y)dy 

JR JR 

u: =? v e ~ z f P(e ?W > u)du = e~ z Ee^ < 00 (2.5) 
Jr >0 

nach (|2.4p . Weiterhin ist jede beschrankte Menge der Form Ct ll ... ) t m (B) £ C in 
einem solchen At tZ enthalten. Daher ist A lokal endlich. 

Satz 2.11. {£(£), t £ M. d } ist Brown- Resnick- stationar genau dann, wenn der 
Poisson-Punktprozess X^ieN^i+^(-) au f C(M. d ) translationsinvariant ist (d.h., 
dass das zugehorige Intensitatsmaft translationsinvariant ist). 
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2.3. Brown- Rcsnick-Stationaritat 



Beweis. Mit den Notationen von oben gilt 

W(v{ti)<yi,...,v(t m )<y m ) 

E f U {/ G c, ( Rd ) : /(*»■) > vj} 1 = 

iSN \j=l,...,m / 

= exp(-A(C tl ,... |tm (B))) (2.6) 
mit B = R m \ Xf =1 (-oc, yj }. 

1st nun der Punktprozess X^Xj+fjM translationsinvariant, so ist A translations- 
invariant, und fur jedes h G M. d folgt nach (|2.6|) 



F(r?(ti) <yi,...,r?(t m ) <y m ) = exp(-A(C tl ,..., fm (B))) 
= exp(-A(C tl+fe) ... itm+/l ( J B))) = P(7j(*i + h)< yi ,..., r){t m + h)< y m ), 

also ist rj(-) stationar. 

Ist andererseits rj(-) stationar, so ergibt sich fur alle h G M. d wiederum nach (|2.6[) 
die Gleichung 

exp(-A(C tl ,... !tm (S))) = ex p(-A(C tl+h ,..., tm+h (B))) (2.7) 
fur alle Mengen 5 = R m \ Xf =l (-oo,yj]. 

Fiir feste ti,...,t m G R d definieren wir nun C tl) ... jtm := {C tl) ... )tro (^4), ^4 G .4)} 
mit .4 := {R m \ X^L 1 (— oo, j/j] : yi, . . . , y m G R}. Da nun „4 stabil unter Vereini- 
gungen ist und aufierdem C tl ,...,t m {Bi) U C tl! ...,t m (B 2 ) = C tl ,...,t m (Bi UB 2 ) fiir alle 
-Sij -B2 G £> m gilt, ergibt sich die Stabilitat von Ct lt ... t t m unter Vereinigung. Mit der 
Siebformel sowie der Gleichung C tl ,...,t m { B i) n C tl ,..,t m (-B 2 ) = C tl ,...,t m {B 1 n J3 2 ) 
folgt nun aus (|2.7p . dass auch der Schnitt zweier Mengen aus Cti,...,tm die Bezie- 
hung 

A(C tl ,..., tm (B)) = A(C tl+h _ tm+h (B)) Vh G R d (2.8) 

erfullt. 

Induktiv so fortfahrend und dabei die Siebformel benutzend erhalt man, dass 
die Gleichung (|2.8p auch fiir jedes Element Ct 1: ... t t m {B) der (nach Konstruktion 
sclmittstabilen) Menge 

M tl ,...,t m ■= {Ci n . . . n C n , n G N, d G Ct u ...,t m } 

gilt. 

Damit gilt die Translationsinvarianz auch fiir die Elemente von 

M:= |J A4 tl ,..., tm 



tl ,...,tn 



und A4 ist schnittstabil. 
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2. Brown- Resnick-Prozesse 



Weiterhin ist A ein Erzeuger von B m und somit erzeugt 

U Ctx,-,^ C X 

meN 

die a- Algebra C. Daher ist M ein schnittstabiler Erzeuger von C. 

Fiir jedes t G R d und B n := R \ (-00, -n] G .4 ist die Folge (Ci(£ n )) ngN C M 
monoton wachsend gegen C(R d ) und wegen (|2.5p gilt 

A(C t (B n )) = A(A t _ n ) < 00 

fiir alle n G N. Damit folgt nach dem Eindeutigkeitssatz iiber Mafie die Transla- 
tionsinvarianz von A auf ganz C. □ 



Satz 2.12. Seien W(-) ein pfadstetiges, intrinsisch stationdres Gaufisches Zu- 
f alls f eld auf R d mit Variogramm 7, Erwartungswert und Varianz <r 2 {-), 
WP ~ u .i.„. W, i G N, j G Z d und T} j) := inf (argsup teRd (w^\t) - ^)) . Es 



2 

gelte 



lim W(t) — = -00 P - /.a. (2.9) 

iii->oo v 2 y 

Weiterhin seien = ^ ^v-OO stochastisch unabhdngige Poisson-Punktprozesse 

i 

auf R (unabhdngig von {W i (•), i G N, j G rmi Intensitat -^e~ x dx, m G N 
sou>ie p G R> - Dann ist 



Z z (t) := max max ( + (t -p-j)- — P ' J; ) , i G 



em Brown-Resnick-Prozess zum Variogramm 7, d. /i. Z3 = Z. 

Dabei ist die Multiplikation p- j als komponentenweise Multiplikation zu verstehen, 
d. h. (pi, . . . ,p d ) • (ji, . . . ,j d ) := (pxjx, . . . ,p d j d ). 

Beweis. Sei ^\t) := W^\t) - 

Zunachst bemerken wir, dass jedes aufgrund von Bedingung (|2.9p P-f.s. end- 
lich ist und auch m\ 3) = sup teRd (Xp ) + $\t)) wohldefiniert ist. 

Wir zeigen nun, dass die Transformation 

$ : C(R d ) -»• R d x C(R d ), X\ j) + eP(-) (T^,X^ + ^(O) 
messbar ist. 

Sei dazu ii,i2>--- e i ne Abzahlung von Q rf . Wir definieren nun fiir G N die 
messbare Abbildung : C(R a! ) — >• R, £ 1— >• Dann ist auch die Abbildung 

M : C(M. d ) — > R, £ h-> sup fc6N Mfc(^) als punktweises Supremum abzahlbar vieler 
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2.3. Brown- Resnick-Stationaritat 



messbarer Funktionen wieder messbar und es ist M(£- ) = . Damit ist wei- 
terhin die Abbildung $ : C(R d ) -»• C(M d ), + ■->■ + - 

messbar. Nun ist nach Konstruktion die „kleinste" Nullstelle von Wir 
definieren fur fc, n £ N die messbaren Abbildungen 

4 n) : C(R d ) -+ R d , rjf (0 = (**' ^ fe) > 

I (oo, . . . , oo), sonst 

und schliefilich r : C(M d ) -> R d , £ ^ lim inf r<; n) (£). 

r ordnet nun jeder (im Unendlichen negativen) Funktion ihre „kleinste" Nullstelle 
zu und ist weiterhin messbar, da punktweise Limites und Infima abzahlbar vieler 
messbarer Funktionen wieder messbar sind. 
Da nun = (to £) gilt, ist <E> messbar. 

Dadurch erhalten wir einen Poisson-Punktprozess 

52 6 (??\x& +&>(■)) 
mit transformiertem IntensitatsmaB 

y(A) = [ \e~ x ¥(x + £ E $- l (A))dx, A£B d xC, 
Jr m 

wobei P das zu W(-) gehorende Wahrscheinlichkeitsmafi sei. Seien nun 

U t : C(R d ) -> C(M d ), f(.)^f(.- t ) 
und 14 : M d x C(R d ) -> M d x C(M d ), (x, /(•)) ^(x + t, /(■ - t)) 

Verschiebungen um ein t G M d . 
Dann gilt 

(* o c/ t )(x^ + #■>(■)) = + - 1)) 



Da £(•) nun nach Bemerkung 12 . 1 Ol Brown- Resnick-st ationax ist, ist das Intensitats- 
maB von fiXi+£i(-) nach Satz 12.91 translationsinvariant (bzgl. Ut) und damit auch 
das IntensitatsmaB von V 5„ci ,(,), , , denn beide Mafie unterscheiden sich nur um 

den Faktor m . Aufgrund der Tatsache, dass <3? mit den Verschiebungsoperatoren 
kommutiert, ist dann auch das transformierte IntensitatsmaB translationsinva- 
riant (bzgl. Vt). 
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2. Brown- Resnick-Prozesse 



Daraus folgt fur jedes j G Z rf : 

max 

c ( mp mp 1 

= m 
te 

^ j) +p-je(p-j 



max 



max 



X?+£\--p-j] 



(2.10) 



Betrachtet man nun jede Seite der Gleichung (|2.10|) fur verschiedene j G Z , so 
erhalt man jeweils stochastisch unabhangige Prozesse. Daraus folgt: 

Z 3 (.) = max max + £{ i} (- - p • j' 



max max 

^ (j) e(p-i-^, P -i+^] 



Weiterhin ist 



max 



max 



(j mod m) i(j mod m) 



'(•) 



(2.11) 



wobei die Modulo-Operation hier komponentenweise zu verstehen ist. Da die ent- 
sprechenden Mengen {i : mod G (p • j — ^f-,p • j + ^r] } disjunkt sind, sind 
die Prozesse auf jeder der beiden Seiten von ()2.1ip fiir ji G Z d , / G N, mit 
jl = j mod m fiir alle I G N jeweils stochastisch unabhangig. Somit gilt auch 



Z^(-) = max max 

T Umo d m), .._™ ip .. + ™| 



A' 



(j mod m) mod m) 



'(■) 



max 

fce{0,...,m-l}° 



max 



modm=fc \T^e(p-j 



max 

iSN 



-Y,P-j+ ,J f] 



max max A > + ( 
fce{o,...m-i} d ieN I 

I I / mp mp 

da [J [p-j-—,p-j + ^- 

j mod m=k 

max max +£?>(■)) =Z(-). 
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2.3. Brown- Resnick-Stationaritat 



Dabei benutzen wir im letzten Schritt, dass Ylke{o m-i} d SieN^x-W a ^ s Summe 
von stochastisch unabhangigen Poisson-Punktprozessen auch ein Poisson-Punkt- 
prozess ist mit Intensitatsmafi X^fce{o m-i} d ~^? e ~ X d ,x = e ~ x dx. □ 

Um andere Darstellungen des Brown-Resnick-Prozesses zu erhalten, verwenden 
wir das nachfolgende Lemma, das eine Folgerung aus Satz 13 aus [KSd09j ist. Der 
Beweis verlauft daher analog. 

Lemma 2.13. SeienW(-) ein pfadstetiges, intrinsisch stationdres Gaufisches Zu- 
f alls f eld auf R d mit Variogramm 7, Erwartungswert und Varianz o~ 2 {-), 

Wi(-) ~ u .i. u . W(-), i G N, sowie Ti := inf (argsup teRd (Wi{t) - ^)), 

Mi := sup teMd (Wi(t) - ^) und schliefilich := Wi(- +T<) - - M t . 

Es gelte (l2~9|) . 



Weiterhin sei durch YlieN^ x i e * n von {^i(')> 2 £ N} unabhangiger Poisson- 
Punktprozess mit Intensitatsmafi e~ x dx gegeben. 

Dann ist X^eN ^ffi^i+Mj,^) ein Poisson-Punktprozess auf R d xRxC(R d ), dessen 
Intensitatsmafi in ein Produkt zerfallt. 

Beweis. Wir definieren zunachst 



£(•) :=W{-) 
M := sup£(t), 



2 ' 



T := inf (argsup (£(i))) 
undF(-) :=£(-+T)-M. 

Um zu zeigen, dass Ylien ^(t^jq+M;,.^) e i n Poisson-Punktprozess ist, beweisen wir, 
dass die Abbildung $ : K x C(R d ) R d x R x C(R d ), (u,Q i-> (T, u + M, F) 
messbar ist. 

Wir wissen bereits, dass die Abbildung £ 1— > T messbar ist (vgl. Beweis Satz 
I2.12|) . Weiterhin betrachten wir eine Abzahlung t\,t2, ■ ■ ■ von Q d . Nun sind die 
Abbildungen /j : C(M d ) — > R, £(•) i-)- i £ N, messbar und somit auch die 

Abbildung auf das punktweise Supremum £ 1— > M = sup igN Ebenso ist auch 
die Abbildung $ : R d x C(R d ) -»• C(M d ), (i, £(•)) ^ £(• + *) messbar, denn fur alle 
ii, . . . , t m , m G N und jede offene Teilmenge -B C M m gilt aufgrund der Stetigkeit, 
dass 

*-\C tl ,... >tm (B)) = R d x P| C tl _ t) ... itfB _ t (S) £ B d x C. 

Somit ist die Abbildung £(•) 1— >■ (T, £(•) — M) 1— )■ F messbar und daher auch 

Damit ist YlieH^{Ti,Xi+Mi,Fi) ern Poisson-Punktprozess auf R d xRx C(R d ) mit 
Intensitatsmafi 

y(A) = [ e~ x ¥(A - (0, x, 0))dx, A £ B d x B x C 
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2. Brown- Resnick-Prozesse 



wobei P das zu (T, M, F) gehorende Wahrscheinlichkeitsmafi sei. 



Nun gilt fur alle y,z G R, i € R , / G C(R ) und ftir jede Borelmenge A C 
R d x R x C(R d ) die Gleichung 



/ e y d9(t,y,f) = [ e~ x [ e*dP(t,y - xj) dx 

JA+(0,z,0) JR JA+(0,z,0) 

= [ e~ x [ e y+z d¥(t, y + z-x,f) dx 
Jr J a 

= [ e~ x ~ z [ e y+z d¥(t,y-xj) dx 
JR J A 



e*dP(t,y-x,f) dx= / eVd*(t,y,f), 
A J A 

also ist e y d9(t, y, f) ein translationsinvariantes MaB in der zweiten Komponente. 
Seien weiterhin 

U t : C(R d ) -+ C(R d ), /(O^ /(•-*) 
und Vt : R d x R x C(R d ) -+ R d x R x C(R d ), (x,y,f(-))^+(x + t,y,f{-)) 

Verschiebungen um ein f G R rf . 
Dann gilt: 

= (Ti + t, Xi + Mi, -t + Ti + t)- Mi) 

= V t (Ti, Xi + M^ -t + Ti + t)- Mi) 

= (V t o*)(Xi +&(.)). 

Da der Poisson-Punktprozess ^ igN 5 CT 2 { .) nach Bemerkung l2.10l und Satz 

12.111 ein translationsinvariantes Intensitatsmafi hat (bzgl. Ut) und die Translati- 
onsoperatoren mit $ kommutieren, ist auch das Bildmafi 9 bzgl. Vt, also in der 
ersten Komponente, translationsinvariant, d. h. 

9(A + (t, 0, 0)) = 9(A) fiir alle A G B d x B x C, t G R d . 

Nun ist 

9([0,l] d x [0,1] x C(R d )) 

= [ e~ x F ( sup (w(t) - ^P-) G [-x, -x + 1], T G [0, l] d j dx 
JR \teR d V 2 / / 

< / e^'P ( sup W(t) > -x) dx = / e x P f sup W{t)>x)dx. (2.12) 

Da sich auBerdem M := sup fg r 0) ud W(t) als Supremum der Gaufischen Zufallsva- 
riablen W(ti), i G N, f, £ (Q fl [0,l]) rf , schreiben lasst und aufgrund der Pfad- 
stetigkeit von W(-) die Beziehung P(|M| < oo) = 1 > gilt, folgt nach Satz 5 



20 



2.3. Brown- Rcsnick-Stationaritat 



aus [LS70| . dass ein e > existiert mit 

E(exp(e|M| 2 )) < oo. (2.13) 

Andererseits ist 

E(exp(e|M| 2 )) = J P (exp(e|M| 2 ) > dz 

= 1 + J P (exp(e|M| 2 ) > dz 

= 1 + J 2ezexp(ez 2 )P (\M\ > z^j dz. 

Es folgt mit (|2.12p die Ungleichungskette 

*([0, l] d x [0, 1] x C(R d )) < f e x ¥ I sup W(t) > x) dx 

\te[o,i] d J 

/— f'CO 
e x dx + I e x ¥(\M\ > x)dx 

< exp + y 2exexp(ex 2 )P (|M| > 

< exp Q + E(exp(e|M| 2 )) < oo nach (l2~T3jl . 

Somit ist das Ma8 e y d^(t,y, /) in den ersten beiden Komponenten translati- 
onsinvariant und jede Menge, deren Projektion auf diese Komponenten im Ein- 
heitswiirfel enthalten ist, hat unter diesem Mafi eine endliche Masse. Also ist fur 
AeC, B £ B d x B durch 



* A {B) := [ eV<W(t,y,f) 
JBxA 



ein translationsinvariantes MaB auf (M. d xl^^x B) definiert, das auf dem Ein- 
heitswiirfel endlich ist. Daher ist dieses MaB ein Vielfaches des Lebesgue-MaBes, 
wobei der (endliche) Vorfaktor von A 6 C abhangt, d. h. ^a(B) = \(A) J B e y dy dt 
mit < X(A) < X(C(R d )) < oo. 
Weiterhin ist durch Q(-) := x ^^d^ e ^ n 

Wahrscheinlichkeitsmafi auf (C(R d ),C) 
gegeben. Mit < A* := X(C(R d )) < oo gilt also 

d9(f,y,f) = X*e~ y dy dt x dQ 

□ 

Satz 2.14. Seien X* und Q wie im Beweis von Lemma \2.13[ Weiterhin sei durch 
II := YlieN^(Si,Ui) e ^ n Poisson-Punktprozess mit Intensitatsmaft X*e~ u du ds ge- 
geben und Fi ~ u .i.v. Q unabhdngig von II. 
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2. Brown- Resnick-Prozesse 



Dann ist 

Z AjX * (t) := max ( U { + F 4 (t - Si) | , t G M d , 

em Brown-Resnick-Prozess zum Variogramm 7 is£, d. h. Z^\*{-) = Z(-). 

Beweis. Seien Xi, Wi, Ti, Mi, Fi wie in Lemma [2. 131 Dann sind nach Lemma [2,13l 
und Lemma fl . 101 sowohl X^ieN ^(t x +m- f ) a ^ s aucn n Poisson-Punktprozesse auf 
R d x R x C(BL d ) mit IntensitatsmaB x e~ u du x Q(dF). Weiterhin lassen sich 
Z(-) und ^4 ; a* (') auf dieselbe Art und Weise aus diesen Punktprozessen gewinnen, 
namlich Z(-) = max !etj T(Ti, Xi + M i: Fi) bzw. Z^\* = maxj eN T(Si, U i: Fi) mit 
r : M. d x R x C(R d ) -> C(R d ), (T, C/,F) ^ [/ + F(- - T). Damit sind Z(-) und 
^4A*(") in Verteilung gleich. □ 



Bemerkung 2.15. Ein ahnliches Resultat erhalten wir, wenn wir die Prozesse 
aus Satz 12.141 auf pL d , p > 0, einschranken. 



Dann ist fur 



7? } = inf (argsup t6pZd (wi(t) - ^ 
M< W = sup ( Wi(t) ° 2{ ^ 
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2 1 \T^\ 



and F i w (-)=Wi(- + 3fO n * - ', t 6 pZ c 



das zufallige MaB y^,- C N^(p) ^ „w ,-,Mn ein Poisson-Punktprozess auf pZ d x 

R x RP zd mit IntensitatsmaB A^V^t x e~"(iu x Q i ~ p \dF) fur ein A( p) > und 
ein Wahrscheinlichkeitsmafi Q^. Eine aquivalente Darstellung Z iX ( p ) von Z er- 
gibt sich analog zu Satz 12.141 Ebenso konnen alle bisherigen Resultate auch auf 
Prozesse mit eingeschranktem Definitionsbereich iibertragen werden. 



Um Z4 zu approximieren, miissen Prozesse mit der Verteilung Q bzw. Q^ p > simu- 
liert werden. Dazu ist es notwendig, diese MaBe explizit zu kennen. Man beachte 
dabei, dass Q im Allgemeinen nicht die Verteilung von W {■ + T) — a 2 {- + T) — M ist 
(und nicht die Verteilung von W(-+T&) - a 2 (- + T^)- M^), Dies erkennt 
man am folgenden Resultat, nach dem die Verteilung von 
ist, falls W(0) = 0. Im Fall einer Brownschen Bewegung W ist aber eine verscho- 
bene Brownsche Bewegung anders verteilt als eine Brownsche Bewegung bedingt 
darauf, dass sich das Maximum im Ursprung befindet. 

Satz 2.16. Sei W(-) wie in Lemma \2. l'J\ und 

T® = inf (argsup tepZd (w{t) - . 
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2.3. Brown- Resnick-Stationaritat 



Wir nehmen weiterhin W(0) = an. Dann ist die Verteilung 

W{-) - | = 0. 

Beweis. Seien A G S(M AZd ) und V G B(R) so, dass 

< / e~ x dx < oo. 
Jv 

Weiterhin sei II = Ylie^^tr^ x +M {p) F (p) 1 ^ er P°i sson "P un ktprozess auf pZ d x 



x MP^ d wie in Bemerkung 12 . 1 51 Dann ist 



q(p\a) = P(n({o} x v x A) = 1 1 n({o} x v x M pZd ) = i), (2.14) 

da das Intensitatsmafi von II ein Produktmafi ist. 

Wenn wir annehmen, dass die Punkte (T^ p \Xi + M- p \ F^) so indiziert sind, 
dass die Folge (Xj)i 6 N realisierungsweise monoton fallend ist (siehe Kapitel 3), 
erhalten wir 

P(n({0} x V x A) = 1 | n({0} x V x R pZd ) = 1) 
= ]T P(T i (p) = 0, Xi + GV\#{i: (if ^X* + M^) G {0} X V} = 1) 

• P(i^ (p) G A I if 5 =0, ^ + Mf } G V, 

#{i : (if \ ^ + m\ p) ) G {0} XV} = 1) (2.15) 

mit 

P(F- (p) G A | if } =0, ^ + Mf 5 G V, 

(f,I 1 +M i (p) )e{0}xF} = l) 

= F(F^ p) G A | if } =0, X, + Aff > G V, 

(T^ p \x ]+ M^)^{0}xV\fj^i) 

= P(lf > G A | if } = 0, X, G V, (if , Xj + uf) i {0} x V Vj / i) 

= P(lf } G A | if } = 0), 

wobei wir ausnutzen, dass Wi stochastisch unabhangig von Xj, Xj und Wj fiir 
alle j i ist. 

Verwenden wir (|2.14p . (|2.15p und 

^ p M?) = Q) X . + M (P) G y | n( { | x y x R pZ d ) = l) = 1, 

so erhalten wir 

qW(A) = p(f. {p) g ^ | Ti = o) = p (w(-) - g a | = 0^ 

fiir alle A G £(M AZd ). □ 
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2. Brown- Resnick-Prozesse 



Beschranken wir uns auf die Annahmen von Lemma I2.13| so konnen wir Q 
als die Verteilung von Fi bedingt auf + M, und beschreiben. Seien II = 
T,ieN 6 (Ti,x l +M l ,F ) i und £ G £(R d x R) so, dass J E e~ x (dt x cfe) G (0,oo). Wei- 
terhin seien N = U(E x C(R d )) und i x ,...,i N so, dass (T ife ,X ifc + M ifc ) G E 1 
fiir A; = 1, . . . , N. Mit Gi, . . . , Gjv bezeichnen wir eine zufailige Permutation von 
i 7 ^, . . . , Fi N , wobei jede Permutation die gleiche Wahrscheinlichkeit habe. 

Satz 2.17. Bedingt auf N = n sind G\, . . . ,G n unabhdngig identisch gemaft Q 
verteilt. 



Beweis. Wir wollen zeigen, dass alle endlich-dimensionalen Rander von G\ , . . . , G n 
Produkte der eindimensionalen Rander mit Verteilung Q sind. Durch Zerlegen 
der Mengen in C und Umindizieren geniigt es, P(Gi G A±, . . . , G ni G A±, G ni +i G 
A2, ■ ■ ■ , G ni+n2 _|_..._|_ n; G A[ I TV = n) equals n!=i Q{Ai) nt fiir paarweise disjunkte 
Mengen A\, . . . , A[ G C, ni, . . . ,ni G N mit n\ + . . . + ri[ < n zu zeigen. Seien 
to = ni + . . . + ri/ und A = Ui=i ^t- Dann gilt 

P(Gi G Ai, . . . , G ni G Ai, G ni+ i G A 2 , . . • , G m G 4j I N = n) 

1 

HGi e Ai,..., G m G Ai I f| n(E x A,) = iV = n) 

ki>ni,...,ki>ni j=l 
k 1 +...+k i <n 

• F(U(E x Aj) = kj, j = 1,... ,1 \ N = n) 
&i fci — n-i + 1 ^2 fcj ~~ n z 

^— ' n n — n — 1 + 1 n — n\ n — m + 1 



ki>ni,...,ki>ni 
ki+...+ki<n 



fcl, . . . , fcj, n - hi - ... - hi) 

■ (1 - Q(A)) n - fcl — ■-*« 

i=l / \fci>ni,...,fc,>n i \ 1 ^ ' < " 1 «/ 

A;i+.-.+fc;<n 

• Q{Ai) kl ~ ni ■ ■ ■ Q(Ai) k '~ n ' ■ (1 - Q(A)) n - kl --~ k '] 

UQ^r) ■ ( E 

\ fei>0,...,fe 

fci+...+fc ; < 

• Q(A x ) kl ■ ■ ■ Q(A) fe ' • (1 - Q(A)) n ~ m ~ kl ~-~ kl 



\i=l / V fei>0,...,fe;>0 

k\+...+ki<n—m 



n — m 
fej, n — m — k\ 



i=l 



□ 
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3 Simulationsmethoden und 
Fehlerabschatzungen 



Aus der Verteilungsgleichheit der Prozesse Z(-), Zi(-), ^(O) ^s(') un d -^4,A*( - ) er- 
halten wir fiinf verschiedene Darstellungen eines Brown-Resnick-Prozesses, aus de- 
nen sich jeweils unterschiedliche Simulationsmethoden ergeben. Dabei tritt stets 
ein Poisson-Punktprozess auf R mit Intensitatsmafi \e~ x dx, A > auf. Daher 
betrachten wir zunachst, wie man diesen simulieren kann. 



3.1 Simulation von Poisson-Punktprozessen 

Sei {N(x),x > 0} ein (homogener) Poisson-Prozess auf R mit Intensitat A, A > 0, 
d. h. N(0) = 0, {N(x),x > 0} hat unabhangige und stationary Zuwachse und 
fur alle x, t > ist N(x + t) — N(x) Poisson-verteilt mit Parameter At (vgl. 
Definition 2.1.1 in |Ros96j ) . Seien Y n , n £ N die zugehorigen „Wartezeiten", d. h. 
Y n := inf{x > : N(x) > n} - inf{x > : N(x) > n - 1}. Dann sind Y n , n £ N 
nach Proposition 2.2.1 aus [Ros96j unabhangig identisch exponentialverteilt mit 
Parameter A. S eicn nun S n '. — 
Dann gilt andererseits 

Lemma 3.1. Es ist II = X^'eN^i e ^ n Poisson-Punktprozess cm/R>o mit Inten- 
sitat A. 

Beweis. Nach Konstruktion gilt fiir < a < b, dass 

n((a,6]) = {#j : Sj G (a, b}} = N(b) - N(a). 

Die Poisson-Punktprozesseigenschaften folgen nun aus der Definition eines Pois- 
sonprozesses: 

1. P(n((o, b]) = n)= P(JV(6) - N(a) = n) = e -A(fe-a) (A(fc-a))" _ 

2. Fiir 44...CK, Ai n A, = 0, i ^ j sind II(Ai), IL(A 2 ), . . . stochastisch 
unabhangig, da {N(x),x > 0} stochastisch unabhangige Zuwachse hat. 

□ 

Die Punkte eines Poisson-Punktprozesses auf R mit Intensitat A lassen sich also 
als Summen von unabhangig identisch exponentialverteilten Zufallsvariablen mit 
Parameter A simulieren. 
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3. Simulationsmethoden und Fehlerabschatzungen 



Aus diesem Lemma erhalt man nun folgendes 

Korollar 3.2. Es ist II = X^jeN ^iog( J-) e * n Poisson-Punktprozess auf R mit 
Intensitatsmafi Xe~' x dx. 



Beweis. Wir stellen zunachst fest, dass die Abbildung 

1 



H : R >0 sh- log 

ein Homoomorphismus ist. Daher ist fur jede Borel-Menge A C R auch H~ 1 (A) 
eine Borel-Menge in R>o und es gilt 11(^4) = Tl(H~ 1 (A)). Auch sind die Urbilder 
disjunkter Mengen unter der Abbildung H disjunkt. Somit bleiben die Eigenschaf- 
ten eines Poisson-Punktprozesses erhalten. 

Zu berechnen bleibt das Intensitatsmafi von IT. Seien a < b G R. Dann gilt 
P(5([a, b}) = 0) = P(n([e- fe , e~ a ]) = 0) 
= exp(-A(e- a - e~ h )) = exp (- J Xe~ x dx\ . 

Also ist das Intensitatsmafi \e~ x dx - wie behauptet. □ 

Einen Poisson-Punktprozess wie in Satz 12.81 konnen wir mit Hilfe des folgenden 
Korollars simulieren, welches ein Analogon zu Lemma 3 aus [Sch02| ist. 

Korollar 3.3. Sei I C R 1 ^ ein endliches Intervall, Ui ~ u .j.u. U(I), i G N. Dann 
ist IT = ^ieN^(iog(— ) u ) e ^ n Poisson-Punktprozess auf R x / mit Intensitatsmafi 

|j| e~ x dx dt. 

Beweis. Folgt direkt aus Korollar 13.21 und Lemma 1 1.101 □ 



3.2 Fehlerabschatzungen 

Im Folgenden werden wir den Brown- Resnick-Prozess gemafi Z(-), Z±(-), ^(O; 
Z^(-) bzw. Z^\*{-) approximieren. Fiir das Berechnen der Approximationsfehler 
(im Sinne der Wahrscheinlichkeit, dass die Approximation nicht der korrekten 
Realisierung des Prozesses entspricht) werden sich die folgenden Lemmata als 
niitzlich erweisen. 

Lemma 3.4. Sei {W(t),t G R} die eindimensionale Standard-Brownsche-Be- 
wegung, Wj(-) ~ u .i.v. W, i G N und YlieN^Xi em davon unabhangiger Pois- 
son-Punktprozess auf R mit Intensitatsmafi Xe~ x dx, A > 0. Weiterhin seien 
tu < < t Q , C,x G R, C > x. Dann gilt 
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3.2. Fehlerabschatzungen 



3i G N : Xi < x, sup Xi + Wi(t) > C 

tdz[t u ,t ] j 

-C 2|^u| mm (\tu\ \ [ i rf, I C — x — \t* 



< Ae-°-^^exp ^ 1-$ 



c-^ V 2 j \ \ ^ 

Beweis. Wir betrachten den Poisson-Punktprozess II := ^ igN ^(Xi,Xi+Wi(0) au ^ 
E x C(M d ) mit Intensitatsmafi 

A(Aj xA 2 ) = [ e~ y F(y + W(-) e A 2 )dy 

fur A\ x A2 £ B xC (vgl. Kapitel 2), wobei P das zu W(-) gehorende Wahrschein- 
lichkeitsmafi sei. 

Fiir f £l definieren wir nun die Menge 

A t * := < (X,W(-)) £Rx C(R d ) : X <x, sup X + W{t) > C \ . 
{ te[o,|f|] J 

Dann gilt: 

Ae" 5 '— / e~ z2 dz dy nach Formel 1.1.1.4 aus [B502] 
-00 V 71 " ./-S=3L 



2\t* 



/2|t*| 



- 00 0FC-y V 2 |i*| / 

00 . _ C y 1 Vgg| / ^ \ , 

Ae e tf — = — exp : — r dy 

c-x y V 2\t*\) 

_ c 2|i*| „/^1\. /"°° 1 „ f V 2 , 2 l**|y I**' 2 



-c 2 i**i ...... (vw r 1 -.._/' fe-i**i) 2 



_ c 2\t*\ f\t*\\ ( /C-x-|t* 



Ae"° — ^ — — exp i- ± 1 - $ 
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3. Simulationsmethoden und Fehlerabschatzungen 



Damit erhalt man die Ungleichung 



P \3i € N : xi < x, sup X + Wj(t) > C 

\ t€[0,\f\] J 

= P (n(^, t1 ) > 0) = 1 - exp (-A(Vl)) 
<A(A| t ,|), dae" x >l-x Vx > 

, c %t*\ (\t*\\ ( „ / C-x-\t*\\\ 



Weiterhin ist 

P ( 3i E N : X < x, sup Xi + Wi(i)>C] 

\ *e[tu,i ] / 

< P ( 3i e N : Xi < x, sup Xi + V^(t) > C j 
V teto,|t«|] / 

+ P ( 3i e N : X < x, sup Xi + iyi(t)>C] 
V *6[0,|t o |] / 

und somit folgt die behauptete Ungleichung. □ 
Dieses Lemma ermoglicht uns eine Aussage iiber das Verhalten von 



P 3f G N : Xi < x, sup Xi + Wi(t) > C 
\ te[o,|t*|] J 

fiir grofie t*. 

Korollar 3.5. Unter den Voraussetzungen von Lemma \3.4\ gibt es L\,L2,to > 
so, dass fiir alle \t*\ > to gilt: 

L i exp I — 

< - log ( 1 - P I 3i € N : Xi < x, sup Xi + Wi(t) > C 

\ \ te[o,\t*\] , 

< L 2 |t*|exp ( ] . 



x 2 / 

Beweis. Seien A, At* wie im Beweis von Lemma [3.41 Dann gilt nach eben diesem: 
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3.2. Fehlerabschatzungen 



Andererseits ist 



r x 9 f°° 2 

A(A t *) = / Xe- y — / e -2 dz dy nach Formel 1.1.1.4 aus [B502] . 
J-oo V 71 " ^-2=1= 

Definieren wir := ip(s)ds mit := exp ( — t) ' * — ^' so S u t 



00 1 a _£ , 
-e 2 

t V2vrs 

1 1 _£ Z" 00 1 s 
-e 2 - 



2vr t 
1 1 



2vr s J 



e 2 



-e 2 - 
2vr t 

> l±-£W) 



1 1 _£ 

e 2 



2vr s J 



+ 

s=t Jt 



1 3 , 

e 2 as 



2tt 
Vt > 0. 



(Eine prazisere Abschatzung, die man durch weitere Iterationsschritte der parti- 
ellen Integration erhalt, findet sich im Beweis von Satz 4 aus [BKS09J.) 

Wegen t~ 3 = o(i~ 1 ) gibt es also ein K\ > so, dass ^>{t) > ^ exp ( — 4 ) fur alle 



t > 1. Da andererseits 



*(*) 



27T<^(t) 



existiert K 2 : = min 4g r 0! i] 



eine stetige Funktion auf M>o und [0, 1] kompakt ist, 
> und somit folgt fiir alle t > die Ungleichung 



Damit erhalten wir fiir \t*\ > 1 die Abschatzung: 



A(A t 



—=eT-dz dy 
nin I K 2 , exp f- (C ^ ] 



2Ae" c ' / e y min A" 



KiJW\ 



C-x 



2\e~°e^r I m in I K 2 , 

lc-x V V 



y 



exp 



2 t* 



y 



dy 



exp 



2|i ! , y 

(y-|*1) 2 
21** I 



f c e^ ' / min ( K 2 \/\ 



Ki\t 
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3. Simulationsmethoden und Fehlerabschatzungen 



. c _JCi /"VFi mrx^ir^ ^ivl**l \ I v 



>2Ae"°e~ / v V|i*|min \ K 2 , — ,— exp -g- 



/| t *i 



-i 



T / t* \ f r , n—r K-At 

> 2Xe~ c exp — min K 2 y/\t^ 



2 J V C-x + |t*|_ 

max{(C - x) 2 , (C - x - y^l) 2 } 



■ exp 



2|t* 



2Ae- c ifi|f| /|t*|\ / max{(C-x) 2 ,(C-x- } 



~ c-x + \t*\ exp v 2 y exp \ 

fur groB genug. 



Mit exp (-%0f. ) "*±?° 1 bzw. exp - ((C -j-f^ )2 exp(-i) und 



P 3z G N : Xi < x, sup X { + Wi(t) > C = 1 - exp(-A(A t .)) 
V *e[o,|f|] / 

folgt die Behauptung. □ 



(i) 

Lemma 3.6. Seien W( ~ W, t £ N, j € N, unabhdngig identisch verteilte 
Standard-Brownsche-Bewegungen auf R. Seien weiterhin X^ieN ^ r ^ 

i 

stochastisch unabhdngige Poisson-Punktprozesse auf R mif Intensitatmaft e~ x dx 
sowie < a\ < a>2 < a>3 < . . . aquidistante Punkte mit Abstand h > und 
di > 4 + /i, oi G /iZ, sowie L > 0, L G /iZ, C G R. Dann gilt: 

-0') _i_ , 1T , / _ 1 \ <- , m/Cj')/ 



3i, j G N : C < Xy> + sup ( W>"(t) - - ) < X\ J> + W> J '(0) 

t£[a,j,aj+L]nhZ 



ie~ c 



— — Z — ex P 



s/a\ — h 



25L [2L I (±VaT=h + y{ai - h) V 25 - log( ^ (ai ~ fe)V25 )) 2 ' 



+ TV7 exp l 2L 



1 / 1 \ / (^7^-2^ 2 

H ;= 1 H ; exp ' 

V2^\ yjai-h-2j 



Beweis. Zunachst definieren wir iir j G N die Poisson-Punktprozesse 



auf C(R rf ) sowie die Mengen 



Aj := I x + w(-) G C{R d ) : C < x + sup ( w(t) - | ) < x + w(0) }> . 
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3.2. Fehlerabschatzungen 



Sei weiterhin P das zu W(-) gehorende Wahrscheinlichkeitsmafi. 

Dann hat II,, j G N, das Intensitatsmafi Aj(^4) := f R e~ x F(W G A — x)dx fur 
A G C (vgl. Kapitel 2). 

Wir wollen nun Aj(Aj) nach oben hin abschatzen. Es gilt: 



C-x < 



sup \W{t) - J ) < VF(0) ] dx 

te[aj,aj+L]nhZ V 2 / / 

[e- x P\C-x< sup ( W(t) - | J < ] dx, 
y te[o^,Oj+L]nhZ V ' ) 



da W(0) = P-f.s. 



< 



e x dx 



+ 



C+ v_2 +log( V_2) 



C 



C - x < sup ( W(t) 

te[aj,a,j+L]nhZ 



< dx 



-c 



exp 



+ e 



-c 



x < sup 

tG[aj,aj+L]nW 



W(t) 



CO] dx 



< e 



-c 



exp 



+ e 



-c 



-log(V-) 



e x P > W(a 



> — 



2 v ■ ? ' 



I W(i) - - ) > x ) dx 

t€[a,j,a.j+L]nhZ V 



+ e 



-log(V0 



e x F{ W{ aj 



< — x fdl 



.r 



I W(t) — - ) > x ) dx 

te[aj,aj+L]nhZ V 



2e~ c 
< exp 



+ e 



+ e 



-c 



-c 



e*F[0>W( ajJ -- 



- > — a/oT ) dx 



e x ¥ W{ aj 



< —\/04 



3 > 



sup W(t) - - > x- + ^ dx, 

te[0,L] V V 2 / 
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3. Simulationsmethoden und Fehlerabschatzungen 



wobei W eine von W unabhangige Brownsche Bewegung sei. 



Nun verwenden wir, dass W(aj) — ~ J\f (— ctj) ist und erhalten damit 



Aj(Aj) < e c — — exp 



+ e 



+ e 



-C 



-c 



-log(V) 



exp 



1 



: exp 



2a.,- 
2a,- 



dx 



sup ( W(i) — - ) > x — y ] dy dx 

te[o,L] 



2e~ c 
< exp 



la i 



+ e 



-c 



e x dx 



exp ( ) dy 



+ e 



-c 



y/2ira. 



■. exp 



(v + £ 



12 

2 



2a,- 



< e 



: exp 



+ e~ C I 1 - exp 



log 



+ e 



-c 



fr- r- * X f 

y-^-iog(^) loo 



2 

v 7 ^ 1 



sup W(t) > x — y dy dx 
.te[o,L] J 



exp [~r ) dy 



i+ 



e 

2 Z" 00 

-= / exp(— z 2 )dz nach Formel 1.1.1.4 aus [BS02j 



: exp 



a j \2 



2a i 



< e 



: exp 



3 \ —C 

1 +e ° 



+ e 



exp 



ex P f ~TT ) d V 



1 ( (2/+t) 2 



2a,- 



exp ( — - — I dy 



yfir(x - y) 



2L 



(3.1) 
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3.2. Fehlerabschatzungen 



Fiir den dritten Summanden aus (|3.ip lasst sich weiterhin abschatzen: 



-c 



-y/*j 1 



exp 



(y + W 

2a,- 



2L 



< e 



Vk(x - y) 
-\f°~j i 



exp 



(x - yf 



2L 



dy dx 



exp 



log(-V-) -/-co 

1 



TT^l- log(^) 



exp 



2a.,- 
2L 



dy dx, 



da x — y > 



log 



2 



log 



< e 



: exp 



2a,- 



• exp 



(g - yf 

2L 



dy dx, 



ist 



da fiir alle x > 1 die Ungleichung x — log(x) > 1 — log(l) = 1 gilt 
und > = 1 ist 



< e 



-c 



3 log(^) 



exp 



fa + f) 2 
2a, 



dy 



Vl 



- log(^ )) 2 
2L 



7T 



exp 



Nun ist die Funktion |x — log(x) auf R minimal an der Stelle x = |. Damit ist 
|x — log(x) > fiir alle x G R, denn es ist || = log(e) > log (|), und diese 
Funktion ist isoton auf R> s . Daraus folgt 



log 



9 3 
-a; + —-x/ci-i 



100 J 5 V 3 
9 



2^ 
10 



> a.,- H — vai — h ( 



100 J 5 

'2 



2 



y/ai-h 5 
2 2 



'3. /a* 





10 


- log ( 






v 2 




- a 



+ 



10 



+ 



log 



2 
log 



log 



/ Vai - a 5 
I 2 2 



log 



Vai ~ 5 
2 2 
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3. Simulationsmethoden und Fehlerabschatzungen 



9 3 , / J(ai - h) V 25 , / - h) V 25 

a,- H — van — ft — log — 

100 J 5 \ 5 & I 2 



K, 



(i) 



+ / y/gj - jj) V 25 _ / y^-ft) V25 



(21 



fiir alle j G N und man erhalt mit K ai := + die Ungleichung 



-c 



i / (y + ^) 2 



^-log(^) 



exp 



2aj 



7T 



cxp 



(^-log(W 
2L 



3 ° [ 1 - exp 



log 



exp 



(1/ + f ) 2 



2a i 



r/y 



-exp -if^ 



• exp 



K, 



a i 



2L 



-r 2L 
< e y — ex P 



iUU 1 • exp 1 



2L 



2L 



Einsetzen in 13.11 liefert schliefilich 
c 2 



Aj(Aj) < e ° exp 



fl.; 



*J \ +e -c 



1+ 



i V2vr 



exp 



dy 



+ e 



-c 



2L 




7T 



exp 



100 

2L 



exp 



2L 



(3.2) 



Nun ist 



3t, j € N : C < + sup ( Wf>(t) - "- ) < + W^O) 
t6[oj,aj+i] 



r{j), 



-0') i w(-J')/ 



1 (3j G N : n, (A,) > 0) = 1 - exp [ - ^ H A j) 



< 



da 1 — e x < x fiir alle x > gilt. 



< 



: exp 



2L 

+ Vv exp 



+ e" 



-c 



loo a i 
2L 



exp 



2L 



nach (13.2 
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3.2. Fehlerabschatzungen 




da der Integrand streng monoton fallend ist 



_ c 8 / y/cTi - h 
<e -exp — 



+ e 



-c 



-C 



< 



/? 



8 exp 



yja\ — h 



+ 



. (8y + 4)^Lexp f -— ) cfy 



200L /2L 

+ — V7 



exp 



A /(ai-/i)V25 



2L 



da y G 



-1 + 



/a V« 
2 ' 2 



a£ [4y 2 ,4(y+l) 2 ] 



< 



exp 



\Ja\ — h 



25L /2L 

+ -9"Vv exp 



- /i + | V( a i - M V 25 - log 



x /(ai-fe)v25 ^ 2 



V2vf V V«i -h-2 



exp 



2L 

(Vol - h - 2) 2 



□ 



Verscharfen wir die Voraussetzungen ein wenig, so erhalten wir 

(i) 

Lemma 3.7. Seien W> ~ W, i € N, j € N, unabhangig identisch verteilte 
Standard-Brownsche-Bewegungen auf R. Seien weiterhin X^ieN ^x (j) ^ r e ^ 

i 

stochastisch unabhangige Poisson-Punktprozesse auf M rmi Intensitatmafi e~ x dx 
sowie < a\ < a>2 < a% < . . . aquidistante Punkte mit Abstand h > und 
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3. Simulationsmethoden und Fehlerabschatzungen 



ai > 16 + h, a± G KL, sowie L > 0, L G KL, C6l. Dann gilt: 



P I 3i, j G N : C < X[ j) + sup ( w{ j) (t) - J J < x\ j) + (0) 

.2e-°( ( , -\ 25L [2L ( (V^h - log(VaT^)) 2 

< — : — exp -y/ai - h + ——a/ — exp ' 



h \ \ v J 9 V 7r \ 2L 

(^7^-4) 2 ' ' 



+ 7^l 1 + Va 1 -/ t -4 )eXP 



Beweis. Der Beweis verlauft analog zum Beweis von 13.61 Da sich jedoch einige 
Konstanten andern und Integrationsgrenzen anders gewahlt werden, wird er hier 
dennoch in Ganze ausgefiihrt. Seien also Uj,Aj,Aj fur j G N genau wie dort 
definiert und P das zu W(-) gehorende WahrscheinlichkeitsmaB. 

Wiederum wollen wir Aj(Aj) nach oben hin abschatzen. Es gilt: 



AjiAj) = f e -*P f ( 



C - x < sup W(t) - - < W(Q) dx 

te[aj,aj+L]nhZ V 2 / / 

C-x< sup ( W(t) - I J < ] dx, 

da W(0) = P-f.s. 

/ e^dx 

'C+^a-+\og{^/o-) 

e -x F l c _ x< gup W ^ _ < dx 

\ te[aj,aj+L]nhZ V ' ) 



C 



-C 



1 



+ e' 



-C 



exp \—\Jaj / 
o 

a7-log(^aj) 



e x P he < sup ( W(t) - - J < j dx 

\ tG[aj,aj+L]nhZ 



< e 



exp ( — -i/aj) + e 



e a; p(0>lV(a j )-^>-2^-, 

aj-log(yaj) 



sup ( W(t) - - ) > x ) dx 

te[aj,aj+L]n/iZ \ 



+ e" 



-c 



f 



e-p ^( ai )-|<-2^aJ, 



sup ( W(t) ~\) > x | dx 
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3.2. Fehlerabschatzungen 



1 ■ 

< exD (-y/aj) + e~ c I e x F (o > W(aj) - ^ > dx 



la ; 



+ e 



aj-\og(^/aj) 



-C 



aj-log(^aj) 



e x P[ W{aj) - ^ < -2 A 



'a , 



3 ' 



te[o,L] 



sup ( - - ) > x- W(aj) + dz, 



wobei W eine von W unabhangige Brownsche Bewegung sei. 



Mit W(aj) - % ~ M {-%,aj) erhalten wir 



hj(Aj) < e c exp (-y/aj) 



+ e 



+ e 



-c 



-c 



■. exp 



2a, 



dy dx 



exp 



P f 

\te[o,L] 



2a,j 



(sup I W(t) — - ) > x — y ] dy dx 
tero.il V 



< exp (— yajj + e 



e x dx 



exp (-^- ) dy 



+ e 



-c 



2^7 1 / (y+f) 2 

exp 1 



2a g - 



/ft; 



1 / 2 

< e~ c ^= exp (^/ajl) + e~ C 



sup W(f) > x — y J dy dx 

1 ( y 2 \, 

exp -— dy 



+ e' 



-c 



- 2 V^I i 



exp 



2a,- 



\/ 27Ta J 

2 f°° 

= / exp(— z 2 )dz dy dx nach Formel 1.1.1.4 aus [BS02] 

7T J x-v 



< e exp + e 



2+ 



o v 



exp [-%r ) ^ 



+ e' 



-c 



^aJ-\og{^/aj) J-oo \J^ a j 

(x - y)' 



- 2 V*J * 1 ( (y + f) 2 

exp 1 



2a,- 



\/n(x - y) 



exp 



2L 



dy 



(3.3) 
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3. Simulationsmethoden und Fehlerabschatzungen 



Fiir den dritten Summanden aus (13.311 lasst sich weiterhin abschatzen 



-c 



a j -\og( y /aJ) J -co ^2ira,j 



i / (y + f) 2 

exp 1 



2a 



exp ( — - — I dy dx 



< e 



-c 



Vk(x - y) 
-V 5 ^ i 



2L 



exp 



aj—logiy/aj) J -co V 3 
1 \[L 



2dj 



IT y/Oj — \Og{y/Oj) 



2L — ' dy dx 



< e 



-C 



Qj-log( v /aJ) J -co 



- 2 V^J e x IL 



7r \ a,- 



exp 



2a,- 



— 2i — ) d y dx > 



< e 



-c 



da fiir alle x > 1 die Ungleichung x — log(x) > 1 — log(l) = 1 gilt 
e x dx I — — exp 



v /a7-iog( v /a7) J -co \J a j 2a j 



— 6XP 2L 



Mit der gleichen Argumentation wie im Beweis von Lemma 13.61 ist die Funktion 
x \— > |x — log(x) positiv und isoton auf M>5 und somit insbesondere auf 



W 0,1 — h, oo) C [4, oo) 



isoton. 



Daraus ergibt sich 



Oj - log (y/oj)) = 

6 



^L> + ^i- log (^) 



25 J 5 V J V 5 
9 6 7 - 



log(^~) + 



2. fa] 



log (Vo) 



2v/(ai -/i) 



log V (ai - ft) 
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3.2. Fehlerabschatzungen 



fur alle j £ N und mit K ai := K$ + ergibt sich die Ungleichung 



-c 



-log(V) 



e x dx I exp 



2a,- 



Va 



7T 



exp 



2L 



< e~ c e° ■ V2tt ■ — exp 

7T 



25 a i 

2L 



exp 



2L 



Einsetzen in 13.31 fiihrt schlieBlich zu 



Aj(^j) < e 



-c 



: exp 



+ e' 



-c 



2L 




TT 



exp 



/a7 + e' 



25 "J 



-C 



2+- 



^ V2vr 



exp 



dy 



2L 



exp 



A" 



2L 



(3.4) 



Nun ist 



P I 3i,j e N : C < + sup I - 7; ) < AT) J; + Wy>(0) 

\ te[aj,aj+L] V 

= P (3j £ N : Uj(Aj) > 0) = 1 - exp I - ^ Aj(A 

\ J6N 

< ^2 Ay(^i)) da 1 - e" x > x fur alle x > gilt 



-0') i 



: exp 



/aj) + 



-c 



2L 

+ Vv exp 



25 a 3 

2L 



exp 



A", 



2L 



nach (EI 



< 



-c 



a\—h 



— =exp + e 



-c 



2+^ V27T 



exp f-^r I 



2L 

+ Vv exp 



<e c — exp (^—\/ai — h^j 



h y„ 2+ v^Ffe 



a\—h 



+ e 



50L /2i 



9/j V vr 



— exp 



_9_ 

25 



(oi - n) 



2L 



exp 



a; 



2L 
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3. Simulationsmethoden und Fehlerabschatzungen 



-c 



< 



2 exp ( — yai — h ) + 



1 



(16y + 16) — exp I — — ) dy 



y 



+ 



50L 2L 




9 V vr 



exp 



{y/ai - h - log(-v/ai - H)) 2 



da y G 



2L 

2 2 



< 



2e 



exp 



/ -\ 25L /2L 

y ai _/ J j + _W_ exp 



ae[V,4( 2 / + 2) 2 
(V«i - ^ - log(Vai - ^)) 2 



2L 



+ 



1 + 



V27r V Va>i -h- 4 



exp 



□ 



3.2.1 Approximation von Z(-) 

Zunachst wollen wir einen Brown-Resnick-Prozess gemafi der Definition 12.41 ap- 
proximieren. Dazu erzeugen wir zunachst unabhangig identisch standardexponen- 
tialverteilte Zufallsvariablen Yi, i£N, und setzen X k := log (^^^—^J. Dann ist 
SieN ^Xi e ^ n Poisson-Punktprozess mit Intensitat e~ x dx nach Korollar 13.21 

Weiterhin sei zu jedem Xj ein (stochastisch unabhangiger) Pfad {Wj(i) — a ^ '■ 
t G R d } gegeben. Definieren wir Z^ k \t) := max i=li ... jfc X { + Wj(t) - fc G 

N, t G M d , so gilt nach Konstruktion zW(t) Z(t) fiir alle t G R d . 

Wir wollen nun den Fehler der Approximation nach k Schritten berechnen. Dazu 
betrachten wir den Spezialfall, dass W(-) die eindimensionale Standard-Brown- 
sche-Bewegung ist. Zudem schranken wir die Approximation auf ein Intervall 
[-6, b] C R, b > ein. Sei nun C := wI teM (ZW '(t) + ^M). 

Proposition 3.8. Seien x(k) < c(k) < 0. Dann gilt 

P(Z ik) (t) ^ Z(t) fiir ein t G [-6,6] | X fc < x(k),C k > c{k)) 

Weiterhin 

P(C fc < c(*)) < 1 - (l - exp (-exp (~^))) ■ (2* (-^) - l)' (3.6) 



und 



\x k > X (k)) < exp( (^jr^ (1 - exp( " ex p(-^ fc ))))- ( 3 - 7 ) 
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3.2. Fehlerabschatzungen 



Die unbedingte Fehlerwahrscheinlichkeit¥(Z^ k \t) ^ Z(t) fiir ein t E [—6,6]) kann 
(Lurch diese Summe der rechten Seiten von (|3.5p . (|3.6[) und ()3.7p beschrdnkt wer- 
den. 

Beweis. Es gilt: 

Z (fe) (i) / Z(i) fiir ein t G [-6,6] | X fc < x(k),C k > c(k) 
<=^ Z ik) {t) < Z{t) fiir ein t E [-6,6] | X k < x(k),C k > c{k) 

^3m>k: X m + W m (t) > Z^ k \t) + > c{k) fur ein t E [-6, 6] 

I X k < x(k),C k > c(k) 

=> 3m > k : sup (X m + W m {t)) > c(k) \ X k < x(k). 
te[-b,b] 

Weiterhin ist die Folge (X m ) m( zf$ nach Konstruktion monoton fallend, d. h. aus 
m > k folgt X m < X k < x(k). 

Daraus erhalten wir folgende Fehlerabschatzung: 



P(Z (fc) (t) / Z(t) fiir ein t E [-6,6] | X k < x(k), C k > c(k)) 

< P(3t E N : X % < x(k), sup {X t + Wi{t)) > c{k)) 

tel-b-b] 

^ , -c(k) b fb\ A ^fC-x-b 

< 4e c(fe) ^T^ — exp - I I 1 - $ 



c{k)-x{k) *\2)\ V \/& 
wobei der letzte Schritt aus Lemma 13 . 41 folgt . 
Weiterhin gilt 



¥(C k > c(k)) > F(Xi > c(fc)/2)P inf W(t) > c(k)/2 

\ie[-6,6] 

= (1 -exp(-exp(-c(A;)/2))) • P f inf^ > c(k)/2 
= (i _ exp( _ exp( _ c(A;)/ 2))) . ( 2 $ - l' 2 



2 



2^6 



nach Formel 1.1.1.4 in jE502] . und 

P(X fc > x(k)) = P(e" Xfe < e _a:(fc) 

x(h) 



t _i exp(-t) 

= ^ P( ^_" 1 1 ) ; 3,W) (l-exp(-exp(-, W ))). 
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3. Simulationsmethoden und Fehlerabschatzungen 



Dann konnen wir die Ungleichung 

P(zW(f) + Z{t) fur ein t G [-b,b]) 
= ¥{Z (k \t) / Z(t) fur ein t G [-6, 6] | X fe < x(k), C k > c{k)) 
■F(X k <x(k),C k >c(k)) 

+ F(Z ( - k \t) / Z(t) ffir ein t G [-6,6] | > x(k)} U < c(k)}) 

■ F({X k > x(k)} U {C k < c(k)}) 

< F(Z {k \t) / Z{t) fiir ein t G [-6,6] | X fc < x(k),C k > c(k)) 
+ F({X k > x(k)} U {C fc < c(fc)}) 

< F(Z {k \t) / Z(t) fiir ein i G [-6,6] | X k < x(k),C k > c(k)) 

+nc k < c(k)) +nx k > x(k)). 

□ 

Um die Konvergenz Z^ gegen Z im Sinne von 

lim F(Z^ k \t) / Z{t) fiir ein t G [-6, 6]) = 

k— >oo 

einzusehen, miissen c{k) und x(k) so gewahlt werden, dass jede der rechten seiten 
von (|3.5|) . (|3.6|) und (|3.7|) fiir k — > oo gegen Null geht. Dies ist z. B. fiir c(A;) = 
-loglog(log(/c)/2) und x(k) = -log(fe)/2 erfullt. 



3.2.2 Approximation von Z 



1 



Nun soli ein Brown-Resnick-Prozess gemafi Korollar 12,71 approximiert werden. 
Wir erzeugen dazu u. i. v. Exp(^)-Zufallsvariablen i G N, j G {l,...,n}, 

und definieren 



4 J) : = lo § 



1 



Dann ist fiir jedes j G {1, . . . , n} nach Korollar 13.21 der Punktprozess X^ieN ^x (3) 

i 

ein Poisson-Punktprozess auf M. mit Intensitatsmafi \eT x ' ■ 

Fiir jedes Paar (i, j) G Nx{l,...,n} sei W i ■ (■) unabhangig identisch verteilt wie 
W(-). 

Seien weiterhin A; G N sowie - wie in Korollar 12. 71 - hi,...,h n G R d gegeben. Dann 
definieren wir 

Z? } (i) := max max + } (i - hj) - ^!£_zM 

j=l,...,n j=l,...,fc 2 

Wiederum gilt lim*^ Z 1 fcl,, " ,fc " ) (t) = Z x (t) fiir jedes t G R d . 

Fiir die Fehlerabschatzung beschranken wir uns erneut auf den Fall, dass W(-) eine 
eindimensionale Standard-Brownsche-Bewegung ist und betrachten die Simulation 
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3.2. Fehlerabschatzungen 



auf einem Intervall [—6, 6] C R, 6 > 0. Weiterhin sei /ti < h>2 < ■ ■ ■ < h n und die 
Menge {hj, 1 < j < n} sei symmetrisch Seien X k := maxj = i v .. in Xu \ C k := 

inf teM (^l*^*) + mi "./ i n " ^J^ ) ■ 

Proposition 3.9. Se« x(fc) < c(/c) < 0. Dann gilt 

F{z[ k \t) ^ Zx(t) fiir ein t G [-6,6] | X k < x(k),C k > c(k)) 

< ic -c(k) b - h c (b-h,)/2 ( 1 _^f c - x -(. b - h iY 



c(k) — x(k) \ \ yjb — h\ 

Weiterhin ist 

¥{C k < c{k)) 
< 1 - (1 - exp(-exp(-c(fc)/2))) 

c(k) - (b-hi) 



V 2^6^67 J V 



l^fb^hx 



und 



F(X k > x{k)) < n fc+l 6Xp( .f V. x{k) \ l - exp(-exp(-nx(k)))). 

(k-iy. 

Beweis. Es gilt 

Z[ k \t) + Z x {t) fiir ein t G [-6,6] | X k < x{k),C k > c(k) 
<=^ z[ k \t) < Zt(t) fur ein t G [—6, 6] | X k < x{k),C k > c(k) 

3j G {1, .. . ,n}, rn > k : + W#)(t - fy) > z{ fc) (t) + > c(fc ) 

fiir ein t G [-6,6] | X k < x(k),C k > c(k) 

=> 3j G {1, . . . , n}, m>k: sup {X$ + W&'^i)) > c(k) \ X k < x(k). 

te[-b-h„,b-hi] 

Nun ist jede Folge (X„ ) ng N monoton fallend und daher folgt aus m > A;, dass 
Xm < Xl — Wir erhalten somit unter Verwendung von —6 — /i n = — (6 — 

hi), dass 

P(zf 1 ""' fcn) (t) / Z x {t) fiir ein t G [-6,6] | X k < x{k),C k > c{k)) 

< ¥(3j G {1, . . . , n}, m > k : sup {X® + Wg>(t)) > c(k) | X k < x(k)) 

te[-b-h n ,b-hi] 

n 

< ^(3? G N : X[ j) < x, sup {x[ j) + wf>\t)) > c(k)) 

j—l te[—b—h„,b—hi] 



< rt^vEIet^ ( 1 + i exp f (c(fc)-^(fc)) 2 - 

V7rc(fc)-x(fe) \ ^ c(k) - x(k) P \ 2(6 -hi) 
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3. Simulationsmethoden und Fehlerabschatzungen 



wobei im letzten Schritt Lemma 13.41 verwendet wurde. 
Fiir c(k) < gilt zudem mit X\ = niaXj = i v „ )T i x\ j \ dass 



\C k > c{k)) > P(X X > c(fc)/2)P inf (W(t) - |*|/2) > c(k)/2 

\te[-(6-/n),6-fti] 



P(Xi > c(fc)/2)P inf (W(t) - 1*1/2) > c(k)/2 
\te[0,b-hi] 



(l-exp(-exp(-c(fc)/2)))- U 



c(k) \/b — h\ 



+ exp(-c(fc)/2) [ 1 - * ( --#= + 



2V& - /ii 

ft) 
2V& - h x 



nach Formel 2. 1.1 A in |BS02| . Weiterhin ist 

P(X^' } > x{k)) = P(e _x * ' < e~ x{k) ) 



o 



fc fc -i exp(-nt) 
* (fc-1)! 



< n * v ( ;_ 1} ; w/ y o m^)* 

= " ^^y" (1 - exp(-exp(-nx(fc)))) 

fiir j = l,...,n. Die Behauptung folgt mit F{X k > x{k)) < *£]=i W{X^ > 
s(ife)). □ 



Konvergenz erhalt man beispielsweise wiederum fiir c(k) = — loglog(log(A;)/2) 
und x(k) = -log(fc)/2. 



3.2.3 Approximation von Z 2 (-) 

Es soil ein Brown- Resnick-Prozess gemafi Kor pilar 12.81 approximiert werden. Wir 
erzeugen dazu unabhangig identisch standardexponentialverteilte Zufallsvariablen 

Yi, i E N, und definieren X k := log ( fc x ) sowie unabhangig identisch auf 
I gleichverteilte Zufallsvariablen Hi, t £ N. Dann ist nach Korollar 13.31 durch 
SieN 3(Xi,Hi) em Poisson-Punktprozess auf R X I mit Intensitatsmafi jyj- e~ x dx dh 
gegeben. Weiterhin seien Wi(-), Wz('), ■ ■ ■ unabhangig identisch verteilt wie W(-). 

Fiir k € N definieren wir 

Zf(i):= max (x % + W t {t - Hi) - ^ ~ ^ 
i=l, ...,fc \ z 
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3.2. Fehlerabschatzungen 



Dann gilt lim^^ zf\t) = Z 2 (t) fur jedes t G R d . 

Fiir die Fehlerabschatzung beschranken wir uns erneut auf den Fall, dass W(-) eine 
eindimensionale Standard-Brownsche-Bewegung ist und betrachten die Simulation 
nur auf einem Intervall [—b, b] C R, b > 0. 

Sei C k = inf te [ a u {Z^\t) + min^g/ a ^ ~ )- Ferner nehmen wir an, dass / sym- 
metrisch sei. SchlieBlich definieren wir t a := b — min^e/ h. 

Proposition 3.10. Unter den genannten Voraussetzungen und Notationen gilt 
P(zf } (i) / Z{t) fiir ein t G [-b,b] | X k < x(k),C k > c(k)) 



< Ae -4k) to exp ( t A( 1 _ q> f C - X - t 



c(k) - x(k) V 2 / V V Vto~ 
Weiterhin ist 

P(C fe < c(k)) 
< 1 - (l-exp(-exp(-c(fc)/2))) 



2y/b-h 1 J V V 1^fb^T 



P(X k > x(k)) < 6XP( /* .^^ (l - exp(- exp(-x(fc)))). 
(fe- 1)! 



Beweis. Es ist 

/ Z 2 (i) fiir ein t G [a, 6] | A fc < C fe > c(fc) 
< ^ 2 (i) fiir ein t G [a, 6] | X k < x(k),C k > c(k) 

==>- 3m > k : X m + W m (t - H m ) > zf\t) + ^ ~ gm) > c(k) 

fiir ein i G [a, 6] | A& < x(k), C k > c{k) 

==> 3m > /c : sup (A m + W m (i)) > c(k) \ X k < x(k) 
te[-t ,t ] 

=> 3i G N : Xi < x{k), sup + Wi(t)) > c(fe). 

te[-t ,t ] 

Die Aussage der Proposition folgt nun direkt aus Lemma 13.41 

Die Abschatzung fiir P(C k < c(k)) lasst sich wie in Proposition 13.91 beweisen. die 
Abschatzung fiir P(X k > x(k)) wie in Propostion 13.81 □ 
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3. Simulationsmethoden und Fehlerabschatzungen 



3.2.4 Approximation von Z 3 (-) 

Wir wollen einen Brown-Resnick-Prozess nach Satz 12.121 approximieren. Wir er- 
zeugen zunachst u. i. v. Exp(^)-Zufallsvariablen i £ N, j £ Z d , und defi- 

nieren 



X® := log 



1 



Dann ist fur jedes j G Z d nach Korollar 13.21 der Punktprozess ^ieN e * n 

i 

Poisson-Punktprozess auf R mit Intensitatsmafi ^e~ x . Fiir jedes Paar G 
N x Z d sei (■) unabhangig identisch verteilt wie W(-). 

Seien weiterhin j min , j moa; G Z d mit j min < j max und k G N gegeben. Dann 
definieren wir 



(t) := _ max _ max XV ; + W> J '{t - pj) 



2 ' 2 J 



if> 6 (- 
Nach Konstruktion gilt 

lim = Z 3 (i) 

fiir jedes t £R d . 

Wir wollen nun auch hier eine Fehlerabschatzung fiir die Approximation durch- 
fiihren. Dabei beschranken wir uns wiederum auf den Fall, dass W(-) eine eindi- 
mensionale Standard-Brownsche-Bewegung ist und betrachten die Simulation auf 
einem Intervall [—b, b] C R, b G pL>o- Weiterhin sei m = 1; j m i n = —jmax S Z<o 
seien so gewahlt, dass 6 + 4 < p • j maa; . 

Zur Vereinfachung betrachten wir den Prozess lediglich auf dem Gitter pL. Wir 
definieren 

T^ p) := inf (argsup ((w^\t) - ^ 
und approximieren Z%{-) durch 

zf ' p) (t) := _ . max . max x\ j) + (t - pj) - ~ pj) 

fiir t G pZ. 

Seien nun X fc := max J=Jmm) ... Jm „ X^, C fc := inf te[ _ 6)6 ] npZ zf ' p \t). 

Unter der Annahme X k < C k gilt fiir jedes j G {j m i n , ■ ■ ■ Jmax}, i > k mit 
f e(-f>f], dass 

max Af > + - f - = X\ j) < X k < C k . 

t£\a,b}n P z 1 2 * 
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3.2. Fehlerabschatzungen 



Somit konnen alle Pfade xjf^ + W^\- — pj) — - - 2 fur i > /c ohne Verschlech 
terung der Approximation abgebrochen werden, sobald X^' < Ck ist. 



Proposition 3.11. Sei c(k) < 0. Dann gilt: 



F(Zf' p '(t) ^ Zf\t) fur ein t £ [-6,6] DpZ |A fc < c(fc),C fc > c(k)) 
16e- c ( fc ) / / vfc " b\ , 100 /&3 



• exp 



(iVpj max -b+ y{pj max -6)V25- io g ( V to — b)v25 -)) 2 
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1 A 1 \ / (VPJmax ~ O - 2) 

1 + / ■ = ~ ex P 



\/2~7T V VPJmax ~ O - 2 

Weiterhin gelten 

nc k < c(k)) 

( ( 2b /p + 1 / c(k)\ f ( p^ 

< 1 - I 1 - exp I exp I — J [1 - exp 

d 4 (i-s>( 



1 - exp(-p/2) \ \ ^86 V 2 

V8b ' V 2 



exp(-c(&)/2) ( 1-$ ( _4=U - ^ 



P(A fc > c(fc)) < (2j max + 1) 6XP( y, C{k) \ l - exp(- exp(-c(fc)))). 

(fc - 1)! 



Beweis. Gemafi obigen Ausfiihrungen haben wir nur den Fehler zu betrachten, 
der durch die Nichtberiicksichtigung der Poisson-Punktprozesse H^' mit j G Z \ 
{jmin, ■ ■ ■ Jmax} entsteht. Wir benotigen also 



in a x 



3(i,j) E N x Z 

c(k) < X\ j) + sup ( W ( J ] (t) - ) < Ap } + (0) 

te[-(pj-a)-(pj-b)}n P z V 2 / 

»(3i £ N, jG Z >inHM) : c(fc) < A^ + sup - < xf 

v te[pj-b, P j+b]npZ \ z J 
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3. Simulationsmethoden und Fehlerabschatzungen 



8e- c <*> ( ( Vpjmaz - l>\ , 100 



exp 



^vS^i |VTW^ - 6) V 25 - log( VM^iZ^))^ 



10 
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1 A , 1 \ ( Wpjmax -b-2f 

H — p= 1 H — = exp 

V27T V VPjmax-0-2J \ 8 

wobei die letzte Ungleichung aus Lemma 13.61 folgt. Analog erhalt man 

F(3(i,j)eNxZ <jmin : 

c(k) < X® + sup (wi j \t) - I) < X{ + Wj(0) 
te[-b-pj,b-pj]npZ \ z / 

. 8e- c <*> / / V-I>-H)min\ , 100 [W 

-— r p l — » — ) + -V7 



• exp 



(^V^^W + |V(-6-pj m i^V25 - log(^ " T l " )V25 )) 2 
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1 A , 1 \ / W-b-PJmin ~ 2) 

1 H — i — exp 1 



\/27T V V~b ~ Pjmin ~ 2 

Die Aussage der Proposition folgt daraus mit 

F(Z^ k ' p \t) ^ Zl p \t) fur ein te[-b,b]D pL \X k < c(k),C k > c(k)) 
< f(3(i,j) G N x Z\ [jniin 5 jmax ] • 

c(fc) < X» + sup (^(t) - |) < Af > + } (0)) . 

Weiterhin definieren wir Zufallsvariablen , i G N durch die Poisson-Punktprozesse 



b/ P 

i£N j=-b/p iei 



Dieser hat das Intensitatsmafi (26/#+l)P(T( p ) = 0)e~ x dx. Da jede Folge (X^) i&i 
monoton fallend ist, konnen wir Y\ > Y2 > . . . annehmen. Damit erhalten wir 

nc k > c{k)) 

> P(Fi > c(k)/2) ■ P ( inf (w(t) -^]>^-\ = 0] 
\te[-2b,2b]n P z \ 2 J 2 1 / 



1 - exp (-(2b/p + I) e - C ( fe )/ 2 P(T^ = 0) 
(l -P Qnf &] (VF(t) - |t|/2) < c(fc)/2^ /P(T^ = 0) 
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3.2. Fehlerabschatzungen 



2&/P + 1 ( c(k)\ / / p 
4 eXp l 2~ J I XP V 2 



l-e-p/ 2 ^ \ v 7 ^ V 2^ 

-exp(-c(ife)/2)(l-*(-^ + -^ 



nach Formel 2.1.1.4 in [B502] . 

Fiir die Verteilung von ergibt sich 

Jmax 



\X k > C(k)) < ]T P (4 > < k )) ^ ^rnax + TO 1J > C(*)). 



J— Jmax 

Daraus lasst sich wie in Proposition ^. 8l die behauptete Abschatzung herleiten. □ 



Setzen wir sogar p • jmin <a—16, b + 16 < p ■ jmax voraus, so erhalten wir analog 
mit Hilfe von Lemma 13.71 die folgende 

Proposition 3.12. Unter den genannten Voraussetzungen und Notationen gilt: 

¥{Zp p) {t) ^ zf(t) fur ein t G [-b, b] n pZ \X k < c(k), C k > c(k)) 
4 e - c ( fc ) 



< — - — I exp (- ^Jpj r , 



100 b 3 ( Wpj max ~b- \og{yJpj max - b)) 

H \ — exp 1 



9 V vr V 46 

Wpjmax - b-Af 



+ ~ /= 1 + / ■ =r 7 ex P 

V27T V VPJmax -0-4/ 

Diese Fehlerabschatzungen weisen fiir b — > oo ein deutlich langsameres Wachstum 
auf als die Abschazungen in den vorhergehenden Abschnitten. Daher konnte es 
fiir grofie Intervalle sinnvoll sein, Z^(-) anstelle von Z(-) zu approximieren. Um 
Konvergenz einzusehen, miissen wir wiederum c(k) und jmax = jmax(k) so wahlen, 
dass jede Fehlerabschatzung aus Proposition 13 . 1 1 1 fiir k — > oo gegen Null geht, z. B. 
c(k) = - \og(k)/2 und j ma x = k. 



3.2.5 Approximation von Z^ x *(-) 

SchlieBlich soil der Brown-Resnick-Prozess auch gemafo Satz 12.141 approximiert 
werden. Sei zunachst I C M d ein Interval! Wir erzeugen nun unabhangig iden- 
tisch mit Parameter A*|/| exponentialverteilte Zufallsvariablen Vi, i € N und 

definieren U k := log ( v . ) sowie unabhangig identisch auf I gleichverteilte 
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Zufallsvariablen Si, i G N. Dann ist nach Korollar 13.31 durch Y2ieN ^(Si,Ui) em 
Poisson-Punktprozess auf Ex / mit IntensitatsmaB X*e~ u du ds gegeben. AuBer- 
dem seien Fi(-), Fz(-), . . . unabhangig identisch verteilt mit Verteilung Q. 

Fur k G N definieren wir nun 

Zf£(t) := .max fc X, + F^t - Si). 

Nach Konstruktion gilt wiederum lim^oo Z±^)(t) = Z&\*(t) fur alle t G R rf . 

Fiir eine Fehlerabschatzung beschranken wir uns erneut auf den Fall, dass W(-) 
die eindimensionale Standard-Brownsche-Bewegung ist und betrachten den Pro- 
zess auf einem Gitter pZ, p > 0. Wir wollen Z 4 A ( p ) auf einem Intervall [—6, 6] 
approximieren. Es sei I =: [v, v] D [—6, b], v G pZ. 

Seien nun X^pn^r v(p)\ ern Poisson-Punktprozess auf ([— v,v] H pZ) x R und 

y^AtzwS, (y)s ein Poisson-Punktprozess auf (pZ \ [—«,«]) x R jeweils mit In- 

tensitatsmafi p\^5 p z(dr)e~ u du und fallend in der zweiten Komponente, d. h. 

Y^ 1 > > . . . sowie y± > y^ > Weiterhin seien F\, fi ~ u .i.v. Q- Dann 

gilt 

Z 4A(P )(-) = max{max(y/ p) + F t (- - ^)), max(yj p) + /,(• - U))) . (3.8) 



Nun definieren wir die diskrete Approximation 

Z[ k £' p \t):= max (y } p) + F t (t - R t 

' i=l, ...,k \ 

fiir t G pZ, fiir welche wir eine Fehlerabschatzung vornehmen werden. 

Weiterhin sei := min ig [__ bib ] npZ Z^^\t). Fiir die folgende Proposition ver- 

wenden wir, dass die Verteilung von W(t) - -| | = ist, wobei VF CHIC 
Standard Brownsche Bewegung und T^> = inf (argsup tgpZ (W(t) — |i|/2)) sei. 

Proposition 3.13. Seien c(k) < und v > b + 4. Dann gilt: 

P(^?(t) / /Sr e^n f G [-6,6] HpZ | < c(fc),C fc > c(*)) 

<64e- c(fc) (l-e-i)- 2 A* 

■(-(-^^(-^-(-^^ 

100v^ / ( 4 v^6 + I y/(v - 6) V 25 - log ^^ 25 )) 5 



H — „ t= exp 
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3.2. Fehlerabschatzungen 



Weiterhin sind 



¥(C k < c(k)) 
< 1 - ( 1 - exp (-X p (2b + p)e~ c(k)/2 



■ ( 1 - (2(o-6)AW) fe - 

^ 1 - e-p/ 2 ^ ^ ^86 V 2 
- exp(-c(fc)/2) 



/ (fc-l)c(fc) N N 

c{k) 



und 



p(y fc (p) > c {k)) < {(2v+ P )\^y 



. exp(-(k-l)c(k)) 
(k-l)\ 



(l_ eX p(- e -(2«+p)AWc W)) _ 



Beweis. Es gilt wegen (|3.8p und X k < c(k) < Ck, dass 



nz[ k ^\t) ± Z^ x(p) (t) fur ein t G [-6,6] HpZ | < c(fc),C* > c(fc)) 
<P[3iGN: sup + W i {t-r i )-\t-r i \/2> c(k) |T (p) = 0] 



< 2 ^ A (p V / e-*p[ 



It - s\ 

sup x + W(t - s) - - — — - > c(k) 
t£[-b,b]npZ 2 



T&> = I dx 



< 2 



P(T(p) = 0) 



( ^ / e" z P ( c(fc) < sup x + W(t) - ^ < x J dx. 
Vseo+pN^ V te[a-M+6]n P z 2 / 



Analog zu Lemma 13.61 ergibt sich daraus 



Hzf'^it) + Z 4xM (t) fur ein i G [-6,6] HpZ \ X k < c(fc),C fc > c(k)) 



4,A(p) 

< 16e" c A*P(T^ = 0)- 1 



exp 



+ 1 + 



a/2tt V yJv-b-2 



exp 



(V^6-2) 2 



100 

90F 



exp 



& + | " b) V 25 - log(^^% 
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AuBerdem ist 
P(T (rt = 0) 

> P^ g( sup ^ (w(t) - M - + ^) < |, - i^i < -|, 

^(p)-¥<-!, sup (w{t)-&-w<p) + &)<!> 



2 ie(p,oo) V 2 2)2^ 

¥(W(p) < 0) • P ( sup W(t) ~l<^ 
\te(o,oo) 2 2 ^ 

2"(l- e », 



nach Formel 2.1.1.4 (1) aus [BSH2] 



1 



4 

Damit folgt die erste Aussage der Proposition. 
Weiterhin gilt 

P(Y fc (p) >x)= P( e - y fe P) < e~ x ) 

\(2v/ P + l) P A^) fc t fc - i eXp( - ( y + | )PA(P) ^ 

(fc-1)! 

< ((2« + p)A ( P ) ) fc exp(-(fc-l)x(A:)) T ( _ (2u + )A (p) t)tft 

[k-iy- Jo 

= ({2v + p)AW) fc eXp( ~ (fc ~ 1):E) (l - exp(- exp(-(2t> + p)\®x))). 

{k-iy. 

Ahnlich wie im Beweis von Proposition 13.111 definieren wir nun Zufallsvariablen 
Y-,Y", ieN, iiber die Poisson-Punktprozesse X)ieN = X^eN *V (p) M s i\< b un< ^ 
£ J6 N <ty = E, £ n ^)l| Si |>6 mit Y{ >Y£>... and If > Y 2 " > . 

i 

Dann gilt 

¥(C k > c(k)) 
>¥(Yl>c(k)/2).¥(Yl' <c(k)/2) 

inf ( Wit) - HP) > °M | TW = 
te[-26,26]npZ V 2 / 2 / 

^1 — exp 

• h - (2(« - 6)A^) fc eXP( (fe _^ } (1 - ex P (- e -(— ,))) 
t kf (W(t) " l*l/ 2 ) < c(fc)/2V/P(TW = 0) J 
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> (l-exp(-A( p )(26 + p)e- c W/ 2 )) 



/ i (fc-l)c(fc) ^ N 

. i _ (2( V - b)X^) k eXH n ± j (l - exp(-e-^- fe ) A(p)c ( fc ))) 



1 - $ 



(k-l)\ 
( c(k) 



1 - e-W2 ^ ^ ^§6 V 2 
-exp(-c(fc)/2) 



□ 



Verschaxfen wir die Voraussetzungen mit 6 + 16 < o, so erhalten wir mit einer 
Rechnung wie im Beweis von Lemma 13.71 die folgende 

Proposition 3.14. Unter den genannten Voraussetzungen und Notationen gilt: 

nz{ k ^\t) + Z^ x(p) (t) filr em t £ [-6,6] HpZ | X k < c(k),C k > c(k)) 
< 16e~ c(fc) (l-e~2)- 2 A (p) 



, A 8 / 1 \ / (V^6-4) 2 

; xpi -^) + 72Tv 1+ 7^iJ exp 

H — _ r- exp 1 



1 V 4& 



Diese Abschatzungen haben fiir 6 — >• oo das gleiche Verhalten wie im vorherigen 
Abschnitt. Daher konnte auch die Approximation von Z&\ besonders fiir grofie 
Intervalle sehr interessant sein. Wahlen wir beispielsweise v = v(k) = A; 1 / 4 und 
c(k) = — log(/c)/4, konvergieren wiederum alle Abschatzungen aus Proposition 
13.131 gegen fiir k — > oo. 



53 



4 Simulationsergebnisse 



Fur eine aktualisierte und ausfiihrliche Simulationsstudie, die insbesondere die 
Aussagen der Satze 12.161 und 12.171 beriicksichtigt sowie die Methoden anhand 
mathematisch scharfere Kriterien vergleicht, sei auf 

Oesting, Marco; Kabluchko, Zakhar ; Schlather, Martin : Simulation Tech- 
niques for Brown Resnick Processes. Submitted to Extremes (2010). 

verwiesen. Ein besonderer Dank gilt dabei den Referees dieses Papers, die mit 
ihren Hinweisen zu einer Verbesserung der Fehlersbschatzungen und der Simula- 
tionsstudie beigetragen haben. 
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